AD-A034  674  PHYSICAL  SCIENCES  INC  WOBURN  MASS 


F/6  20/5 


ANALYTIC  MODELING  OF  ELECTRICALLY  EXCITED  D2/HCL  AND  HCL  LASER  — ETC(U) 


unclassified 


JUL  76  R TAYLOR*  G CALEDONIA*  P LEWIS*  P 
PSI-TR-58 


N00014-75-C-0035 

NL 


1.0 

!:  him 
m 

HIM 

|2! 

1 11 

t m 
t ia 

u 

bV.U 

1= 

2 ~ 

Ilia 

m 

1.25 

JA 

MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARDS-1963-j$ 

I 

i 


R.  TAYLOR,  G.  CALEDONIA,  P.  LEWIS,  P.  WU, 
J.  D.  TEARE  AND  J.  CRONIN 


Sponsored  by 

OFFICE  OF  NAVAL  RESEARCH 
Department  of  the  Navy 
Arlington,  Virginia  22217 

Contract  No.  N00014-75-C-0035 
NRL  Req.  173-75-RQ-06009 

1 

! 


D D C 

iprarPMUE 


JMI  81  19TT 

tEcsEtnns 


DISTRIBUTION  _ STATEMENT  A 

Approved  for  public  release; 
Distribution  Unlimited 


ANALYTIC  MODELING  OF  ELECTRICALLY  EXCITED 
D /HC1  AND  HC1  LASER  EXPERIMENTS 


FINAL  REPORT 


by 


R.  Taylor,  G.  Caledonia,  P.  Lewis,  P.  Wu, 
J.  D.  Tea  re*  and  J.  Cronin 


July  1976 


Sponsored  by 

OFFICE  OF  NAVAL  RESEARCH 
Department  of  the  Navy 
Arlington,  Virginia  22217 

Contract  No.  N00014-75-C -0035 

NRL  Req.  1 73 -75 -RQ -06009 


PHYSICAL  SCIENCES  INC. 
30  Commerce  Way 
Woburn,  Mass.  01801 


^Present  address:  W.  J.  Schafer  Associates , Inc.  , Wakefield,  MA  01880 

Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government 


ABSTRACT 


’1 


. 

An  analytic  capability  has  been  developed  to  model  two  type  of 
laser  experiments:  (i)  a D2/HC1  EDL  and  (ii)  an  e-beam  sustainer 
discharge  HC1  laser. 

For  the  former  device,  a discharge  in  D2/inert  gas  mixtures 
provides  excitation  of  the  vibrational  mode  of  D2.  This  gas  is  then 
supersonically  expanded  and  HCi  injected.  The  gases  mix  and  flow 
into  a cavity  where  vibrational  energy  transfer  pumps  the  HCl.  Com- 
plete modeling  of  this  sequence  of  steps  has  been  achieved.  The 
relevant  electron  impact  rate  constants  for  D2  (and  f^)  as  well  as 
the  vibrational  energy  transfer  rate  constants  for  the  D2/HC1  molecular 
system  were  selected  from  available  data  and  theory.  Comparison  of 
the  results  of  this  model  with  data  from  the  NRL  EDL  experiment  were 
in  qualitative  agreement,  although  more  definitive  comparison  was 
precluded  by  experimental  difficulties.  A conceptual  design  of  a next 
generation  experiment  was  performed  and  indicates  the  requirement  for 
an  e-beam  sustainer  discharge  operating  at  about  one  atmosphere  plenum 
pressure.  Under  these  conditions  laser  demonstration  of  a D2/HCl  EDL 
is  predicted  to  be  feasible. 


In  the  second  experiment,  HCl/inert  gas  mixtures  are  directly 
excited  by  a pulsed,  high  current  e-beam  sustainer  discharge.  The 
discharge  conditions  were  modeled  by  developing  a set  of  electron  impact 
cross  sections  for  HCl  by  a critical  review  of  existing  data.  The  results 
of  model  calculations  indicate  that  for  the  conditions  of  the  NRL  experi- 
ment, sufficient  atoms  may  be  generated  in  the  discharge  to  prevent 
attainment  of  high  gain.  Definitive  modeling  of  the  atom  generation  is  not 
possible  at  present  due  to  lack  of  kinetic  information.  A significant  range 
of  experimental  variation  could  not  be  achieved  due  to  device  limitations. 
Additional,  controlled  experiments  are  suggested  to  assess  the  critical 
technical  questions  concerning  an  HCl  discharge  laser. 
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I.  D2/HC1  EDL  MODEL 


A.  INTRODUCTION 

This  Section  of  the  report  describes  analytical  modeling  of  the 
D2/HCI  electric  discharge  gas  dynamic  laser  (EDGDL)  under  investiga- 
tion at  NRL.  The  concept  of  this  device  is  as  follows:  D2/diluent  mixtures 
are  passed  through  a low  to  medium  pressure  D.  C.  glow  discharge. 

In  principle,  a high,  i.  e.  , atmospheric  pressure  sustained  discharge 
could  also  be  employed.  It  is  expected  that  at  the  proper  value  of  E/N, 
a substantial  fraction  of  the  discharge  energy  will  be  dissipated  in  exciting 
the  vibrational  mode  of  the  D2.  The  vibrationally  excited  D2  is  then 
expanded  through  a supersonic  nozzle  and  mixed  with  HCl.  Further  down- 
stream in  the  cavity,  the  D2  transfers  its  vibrational  energy  to  the  HCL 
via  a V-*V  process: 

D2  (v  = 1)  + HCl  (v  = 0)  -»  D2  (v  = 0)  + HCl  (v  = 1)  (I-  1) 

In  addition,  rapid  intramode  V-»V  pumping  occurs  among  the  vibrationally 
excited  HCl  molecules: 

HCl  (v  = n)  + HCl  (v  = 1)  ->  HCl  (v  = n + 1 ) + HCl  (0)  (1-2) 

The  occurrence  of  processes  (I - 1 ) and  (1-2)  under  conditions  of  high  vibra- 
tional and  low  translational  temperatures  is  an  example  of  Teare -T ra inor  , 
pumping  and,  in  principle,  can  produce  a partial  inversion  on  the  HCl  upper 
vibrational  levels. 

A schematic  of  the  NRL  experimental  apparatus  is  shown  in 
Fig.  1-1.  The  discharge  is  struck  between  a set  of  ballasted  pins  and  the 


Further  details  concerning  the  NRL  EDGDL  and  the  results  of  the  experi- 
ments performed  on  this  apparatus  may  be  found  in  the  following  series  of 
reports:  ARPA-NRL  Laser  Program  Semi-Annual  Technical  Reports: 

1 July  1973-31  Dec.  1973  NRL  Memorandum  Report  2846,  July  1974  p.  19; 

1 Jan.  1974-30  June  1974,  NRL  Memorandum  Report  3005,  April  1 975, 
pp.  52-56;  1 July  1974-31  Dec.  1974,  NRL  Memorandum  Report  3084, 

July  1975,  pp.  18-27;  1 Jan.  1975-30  June  1975,  NRL  Memorandum  Report 

3217,  Feb.  1973,  pp.  13-21. 


EPOXY  SCREEN  SUPPORTS  PINS 
SERVES  AS  TURBULENCE  GENERATOR 


end  of  the  nozzle  blades.  The  discharge  length  is  about  20  cm.  Typical 
plenum  conditions  are  100  - 150  torr  of  30%  D2  - 70%  He  mixture.'  The 
discharged  D2  is  expanded  to  approximately  Mach  2 and  mixed  with  a 
small  amount  of  HC1.  Final  gas  mixture  composition  in  the  cavity  is  typically 
0 3%  HCl/30%  D2  /70%  He  at  a pressure  of  ~6  torr.  Observation  windows 
extend  from  about  the  end  of  the  nozzles  to  15  cm  downstream.  An  addi- 
tional 14  cm  of  downstream  observation  is  available  by  repositioning  the 
windows.  The  transverse  cavity  dimension  is  25  cm. 

The  primary  diagnostics  consist  of  the  discharge  parameters,  gas 
flow  rates  and  pressures,  and  HC1  infrared  fluorescence.  This  latter 
measurement  is  made  via  an  optical  system  which  translates  parallel  to 
the  cavity  during  a run.  Data  of  HC1  fluorescence  intensity  versus  down- 
stream distance  are  obtained. 

In  the  following  sections  of  this  report  we  will  discuss  the  analytical 
modeling  of  this  experiment,  comparisons  with  data,  and  preliminary 
concepts  for  an  improved  device.  Because  of  the  geometry  of  the  experi- 
ment, the  modeling  naturally  divides  into  major  tasks:  (i)  the  discharge 

or  plenum  region  and  (ii)  the  cavity  region.  A major  portion  of  the  first 
effort  was  devoted  to  the  development  of  a state-of-art  electron  impact 
cross  section  data  package  for  D2  over  the  electron  energy  range  of  interest 
in  this  study.  Because  of  the  existing  limited  data  base  for  D2,  cross 
sections  were  first  developed  for  H2,  and  theoretical  arguments  used  to 
extrapolate  the  H2  data  to  D2 . Under  the  cavity  modeling,  we  will  discuss 
the  vibrational  energy  transfer  rate  constants  for  the  D2/HCl  molecular 
system,  approximate  mixing  models,  and  calculations  of  fluorescence  an 
small  signal  gain  under  typical  experimental  conditions.  Finally,  some 
limited  comparison  with  data  are  discussed.  It  is  shown  that  fluid 
mechanical  problems  associated  with  the  device  during  much  of  the 
program  preclude  a detailed  comparison.  Suggestions  are  developed  to 
construct  a second  generation  experiment  by  which  to  provide  laser  device 

verification. 


* More  specific  experimental  conditions  are  presented  later  in  this 
report. 
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B.  DISCHARGE  MODELING  IN  H2/D2 

1.  Boltzmann  Code 

The  purpose  of  this  task  was  to  determine  the  electron  energy 
distribution  and  molecular  excitation  rates  arising  from  electric 
discharges  in  the  gases  Hz  and  D2  with  the  ultimate  goal  of  providing 
laser  performances  predictions.  The  phenomena  modeled  involve:  (i)  the 
application  of  an  electric  field  across  a partially  ionized  gas  which  causes 
the  electron  acceleration;  (ii)  the  various  inelastic  collisions  of  electrons 
with  the  gas  molecules;  (iii)  the  resulting  electron  energy  distribution 
caused  by  the  combination  of  these  inelastic  collisions  with  the  electric 
field  strength;  and  (iv)  the  excitation  rate  constants  determined  by 
integrating  the  relevant  excitation  cross  sections  over  the  resulting 
electron  velocity  distribution. 

The  electron  energy  distributions  in  discharges  are  generally 
non- Maxwellian  and  must  be  determined  by  a solution  of  the  Boltzmann 
transport  equation.  Fortunately,  much  effort  has  been  directed  toward 
the  solution  of  this  problem.  The  computer  code  used  in  the  present 
analysis  was  originally  developed  by  Carleton  and  McGill*  and  has  since 
been  updated  by  Lenander^  among  others.  The  inelastic  processes 
considered  and  included  in  the  code  are:  (a)  momentum  transfer,  (b) 
rotational  excitation,  (c)  vibrational  excitation,  (d)  dissociative  attachment, 
(e)  molecular  dissociation,  (f)  electron  excitation,  and  (g)  ionization.  The 
main  thrust  of  the  present  effort  was  to  develop  a cross  section  data 
package  for  the  H£/D£  molecular  system  to  be  used  with  the  existing 
computer  code. 

The  derivation  of  the  Boltzmann  equation  is  discussed  in  detail  in 
Carleton  and  McGill.  * The  specific  solution  starts  with  the  Boltzmann 
equation  for  a spatially  uniform  electron  distribution  f(v,  t): 

§ + i • ’ (f ) i'-3’ 

where  a is  the  acceleration  due  to  the  driving  force,  V is  the  gradient  in 
velocity  (v)  space,  and  t is  time.  A solution  is  assumed  of  the  form: 

f W = fQ  (x)  + £ * 1 (v)  + x £)  • v f2  (v)  (1-4) 
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where  and  .B  are  the  electric  and  magnetic  fields,  respectively. 

An  expansion  of  f (v)  in  spherical  harmonics  in  v space  results  in 
the  following 
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E is  the  rms  amplitude  of  the  electric  field  with  angular  frequency  uu  ; 
g (u)  = f (u);  u (eV)  = kv^  where  k is  a unit  conversion  factor;  e/m  is  the 
electron°charge  to  mass  ratio;  Nj  is  the  density  of  species  i with  mass 
which  has  elastic  cross  section  0^c  (u),  an  average  rotational  cross 
section  Oir,  and  inelastic  cross  section  o^j  for  inelastic  process  j with 
threshold  energy  Ujj  (eV).  Wj  is  the  average  energy  loss  in  a rotational 
collision  with  species  i. 

In  Eq.  (1-5)  it  can  be  seen  that  the  first  term  represents  the 
momentum  transfer,  the  second  term  the  rotational  excitation  collisions, 
and  the  third  term  includes  all  other  inelastic  processes  such  as  vibrational 
excitation,  ionization,  dissociation,  electronic  excitation,  and  dissociative 
attachment.  In  our  calculations  the  rotational  cross  section,  CJjr,  is  set  to 
zero  and  rotation  is  treated  as  another  group  of  inelastic  processes. 

The  expressionfor  f^  (u)  and  f^  (u)  are  given  in  terms  of  g (u)  as: 

+ - ^ ■““>  £ ”-61 

f2  <">  = £1  <">  (I-7) 
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where  0 = (e  B/m)  and 
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The  small  perturbation  solution  as  described  above  is  valid  for  a 
periodically  fluctuating  electric  field  _E,  with  frequency  uu,  and  magnetic 
field  J3.  It  is  of  questionable  accuracy  at  high  field  strengths  and/or  low 
frequencies;  however  it  is  valid  for  uu  = 0.  For  the  current  application  con- 
siderable simplification  exists  because  the  electric  field  is  constant,  there- 
fore uu  = 0,  and  there  is  no  magnetic  field.  Hence,  Eq.  (1-8)  becomes 
a (u)  = 1/v,  f2  (u)  = 0 and  further  simplification  of  the  other  relations  are 
obvious.  Further  discussion  of  the  limits  of  validity  of  the  theory  may  be 
found  in  Ref.  2. 

2.  Data  Base 


The  only  comparison  between  predicted  (i.  e.  solutions  developed 
from  Eq.  (1-3))  and  measured  electron  transport  properties  in  H2/D2 
discharges,  which  has  been  performed  for  a wide  range  of  the  Townsend 
parameter  E/N,  is  the  study  of  Englehardt  and  Phelps.  ^ At  the  time  of 
this  latter  study  the  available  cross  section  data  base  for  elastic  and 
inelastic  collisions  between  electrons  and  H2/D2  was  very  limited  and 
many  of  the  cross  sections  were  deduced  by  Englehardt  and  Phelps 
through  a technique  involving  trial  and  error  comparisons  between  predicted 
and  measured  transport  properties.  (In  this  analysis  theoretical  and 
experimental  guidelines  were  used  wherever  possible  in  specifying  the 
cross  sections. ) Since  the  time  of  this  analysis  a considerable  amount  of 
new  data  concerning  the  collision  processes  in  the  system  have 

been  amassed,  and  it  has  become  evident  that  the  Englehardt  and  Phelps 
cross  section  base  requires  substantial  revision. 

This  re-interpretation  is  due,  in  part,  to  the  fact  that  the  measured 
transport  properties  of  H2/D2  are  now  more  accurately  defined.  In 
particular,  it  has  been  found  that  the  transport  data  used  in  Ref.  3 for  H2 
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at  large  E/N,  > 6 x 10  V/cm  , were  erroneous.  This  led  to  an 
incorrect  definition  of  the  momentum  transfer  and  electronic  state  excita- 
tion cross  sections  at  high  electron  energies,  i.  e. , > 10  eV.  More 
important,  from  the  viewpoint  of  laser  performance  predictions,  is  the 
fact  that  the  theoretical  treatment  applied  by  Englehardt  and  Phelps3  to 
specify  the  cross  sections  for  rotational  excitation  has  since  been  found  to 
be  inadequate,  and  therefore,  both  the  rotational  and  vibrational  cross 
sections  presented  in  Ref.  3 are  invalid.  In  the  present  analysis  the  cross 
section  data  base  has  been  updated.  The  goal  of  this  update  was  not  to 
provide  the  most  "exact"  comparison  with  the  available  transport  data,  but 
rather  to  ensure  that  the  comparison  between  predicted  and  measured 
transport  properties  was  sufficiently  good  so  that  laser  performance 
predictions  could  be  developed  with  a reasonable  degree  of  confidence. 

Much  of  the  relevant  cross  section  data  base  for  H 2 molecules  may 
be  found  in  the  recent  data  compilation  by  Kieffer;4  however,  no  recom- 
mended cross  sections  nor  comparison  with  theoretical  predictions  are 
included  in  this  compilation.  Therefore,  a search  of  the  recent  literature 
was  performed,  the  available  experimental  data  and  theoretical  predictions 
for  each  process  were  reviewed,  and  a "best"  cross  section  representation 
chosen.  These  selected  cross  sections  represent  the  authors'  personal 
judgement  (and  perhaps  biases)  concerning  the  relative  merits  of  the 
various  measurements  and  theoretical  calculations.  In  most  cases  the 
high  energy  (>  10  eV)  portion  of  the  cross  sections  had  to  be  extrapolated 
as  no  theoretical  predictions  or  measurements  are  available. 

The  most  studied  inelastic  process  for  H2  is  the  cross  section  for 
the  v = 0 1 vibrational  excitation.  The  data  base3*  for  this  process 

is  shown  in  Fig.  1-2.  The  larger  cross  sections  were  deduced  from  drift 
tube  results,  while  the  remainder  of  the  data  were  determined  by  electron 
beam  experiments.  The  latter  measurements  are  preferred,  since  the 
data  analysis  for  a beam  experiment  is  more  straightforward  than  for  the 
drift  tube.  The  cross  section  selected  to  be  used  for  the  present  calcula- 
tions was  essentially  that  of  Ehrhardt  et  al*9  extrapolated  into  the  higher 
energy  data  of  Trajmar  et  al.  9 The  data  for  excitation  into  higher  vibra- 
tional levels  of  H2  is  more  sparse,  and  the  choice  of  these  cross  sections 
is,  therefore,  considerably  more  subjective.  As  an  example,  the  data  base 
for  vibrational  excitation  to  v = 2 is  shown  in  Fig.  1-3.  Again  the  cross 
section  selected  was  that  of  Ehrhardt  et  al^  extrapolated  into  the  high 
energy  data  of  Trajmar  et  al.  9 Fortunately,  the  larger  excitation  cross 
section  is  to  the  first  vibrational  level. 
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Rotational  excitation  was  handled  in  an  approximate  manner  in  the 
present  calculations.  Both  H2  and  D2  are  composed  of  ortho  and  para 
states  with  certain  statistical  weights.  The  relative  concentrations  of 
each  rotational  state  are  also  a function  of  temperature.  Cross  section 
data  for  the  J = 0 ■*  2 and  1 -*  3 transitions  are  well  defined  for  both  H2  and 
D2  except  at  the  higher  electron  energies  in  D2.  11  Little  information  is 
available  on  higher  transitions  such  as  J = 2-*  4.  The  existing  version  of 
the  Boltzmann  code  uses  a continuum  rotational  excitation  model,  but  this 
is  not  appropriate  for  molecules  such  as  H2/D2  which  have  large  rotational 
spacings.  Therefore,  the  code  has  been  modified  to  calculate  each  rota- 
tional transition  as  a separate  inelastic  process.  The  cross  section  for 
each  rotational  transition  must  be  put  into  the  code  separately,  properly 
weighted  by  the  population  of  that  rotational  state.  In  the  present  predic- 
tions the  H2  and  D2  molecules  were  taken  to  be  50%  in  the  J = 0 state  and 
50%  in  the  J = I state.  As  discussed  in  Appendix  I-A,  this  turns  out  to  be 
a very  reasonable  approximation. 

For  the  remaining  cross  sections  for  H2  the  following  correlations 
were  used:  (a)  The  momentum  transfer  cross  section  deduced  by  Gibson11 

was  applied  at  low  electron  energies,  (b)  The  dissociation  cross  section 
used  was  that  determined  by  Corrigan,  ^ albeit  corrected  by  more  recent 
measurements  of  the  ionization  cross  section.  (3)  Lastly,  the  high  energy 
momentum  transfer  and  electronic  state  excitation  cross  sections  deter- 
mined by  Englehardt  and  Phelps3  were  used.  Although  it  is  recognized 
that  the  latter  cross  sections  are  incorrect,  these  processes  are  only 
important  at  characteristic  energies  beyond  those  for  optimum  vibrational 
excitation.  Since  such  discharge  regimes  are  not  of  interest  for  the  present 
laser  applications,  it  was  felt  that  the  effort  required  to  develop  accurate 
cross  sections  in  this  energy  region  was  not  warranted. 

The  data  base  for  inelastic  collision  cross  sections  for  D2  is  very 
limited,  however,  it  is  expected  to  be  analogous  to  that  for  H.,.  3 Where 
data  were  not  available  for  D2,  cross  sections  were  constructed  which 
were  similar  in  shape  and  magnitude  to  those  for  H2  but  adjusted  to  have 
energy  onsets  appropriate  to  D2.  The  onset  behavior  of  the  v = 0 -»  1 
vibrational  excitation  cross  section  for  D2  has  been  specified  by  Gibson.  11 

The  appropriate  cross  sections  of  interest  for  H2  are  shown  in 
Fig.  1-4.  The  data  for  D2  are  quite  similar.  For  simplicity,  only  one  of 
the  cross  sections  for  rotational  and  vibrational  excitation  is  shown.  It  is 
to  be  remembered  that  each  of  the  processes  in  Fig.  1-4  results  in  a 
distinct  energy  loss  to  the  electrons  and,  thus,  will  affect  the  electron 
distribution  function  in  a different  manner. 
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Discharge  predictions  have  been  made  for  both  pure  H2  and  D2  over 
the  E/N  range  of  10’16  - 10“15  volts-cm2.  As  a check  on  the  validity  of 
the  results,  comparisons  were  made  between  the  predicted  electron  drift 
velocity  and  characteristic  electron  energy  and  the  best  available  drift  tube 
measurements  of  these  quantities  as  presented  in  Ref.  13.  The  comparison 
for  H2  is  shown  in  Fig.  1-5.  Although  it  is  obvious  that  further  refinements 
could  have  been  made,  the  agreement  is  quite  reasonable,  with  differences 
between  predictions  and  data  being  7%  at  the  worst  for  E/N  between  1 and 
6 x 10~*°  V-cm2.  The  difference  becomes  worse  at  increased  E/N  as  a 
result  of  the  improperly  specified  cross  sections  at  higher  electron  energies. 
In  particular,  the  predicted  characteristic  energies  are  too  high  and  drift 
velocities  too  low,  continuing  the  trend  which  can  be  observed  in  Fig.  1-5. 
Conversely,  the  predictions  at  higher  E/N  are  in  fair  agreement  with  those 
developed  by  Enhglehardt  and  Phelps.  3 This  suggests  that  the  changes  made 
in  the  low  energy  cross  sections,  i.  e.  rotation,  vibration,  do  not  affect  the 
predicted  transport  properties  at  high  E/N. 

Such  a comparison  cannot  be  made  in  the  case  of  D£  inasmuch  as 
there  are  no  measurements  of  the  room  temperature  transport  properties*3 
for  E/N  ^ 10"*^  V-cm2.  Gibson**  has  shown  that  the  low  energy  cross 

sections  for  D2,  which  are  essentially  identical  to  those  used  in  the  present 
analysis,  provide  good  agreement  ( < 1%)  with  the  observed  transport 
properties  of  D2  in  the  range  of  E/N  of  1 to  10  x 10" *7  V-cm2.  Nonetheless, 
he  points  out  that  the  rotational  cross  sections  cannot  be  determined 
uniquely  from  the  available  electron  swarm  data.  Inasmuch  as  the  conjec- 
ture that  the  H2  and  D2  cross  sections  are  similar  cannot  be  confirmed  at 
present,  the  predictions  for  D2  presented  in  the  following  pages  must  be 
considered  as  qualitative  rather  than  quantitative. 

Predictions  of  the  energy  deposition  and  excitation  rate  constants  for 
pure  H2  are  shown  in  Figs.  1-6  and  1-7  respectively.  While  the  excitation 
rate  constants  are  required  to  perform  kinetic  calculations,  the  relative 
energy  distribution  can  be  used  to  determine  optimum  discharge  parameters. 
As  can  be  seen  from  Fig.  1-6,  optimum  pumping  into  the  vibrational  mode 
of  H2  occurs  at  E/N  ~ 3x  10"*^  volts-cm^.  At  lower  E/N  rotational 
excitation  begins  to  dominate,  while  at  higher  E/N  dissociation  becomes 
important.  Note  vhat  even  at  the  optimum  operating  condition  only  ~ 70% 
of  the  electron  enc  rgy  loss  goes  into  vibration  with  the  remainder  going 
into  translational  and  rotational  heating.  From  Fig.  1-7  it  can  be  seen  that 
the  vibrational  excitation  rate  constants  are  also  near  their  maximum  at 
an  E/N  of  3 x 10“*^  volts  -cm^  and  decrease  considerably  at  lower  E/N. 

Note  that  the  excitation  of  the  first  vibrational  level  of  H2  is  by  far  the 
dominant  process. 
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Fig.  1-5 


Comparison  Between  Predicted  and  Measured  Discharge 
Properties  of  Pure  H2  at  T = 293°K. 


Fig.  1-6  Calculated  Percent  Discharge  Energy  Deposition  int 

Various  Internal  Modes  of  Pure  H,. 


From  these  calculations  it  was  found  that  the  dissociative  attach- 
ment cross  section  was  sufficiently  small  so  that  this  process  played  no 
important  role  in  the  discharge.  Furthermore,  over  the  range  of  E/N 
considered,  the  electronic  excitation  and  ionization  processes,  although 
instrumental  in  defining  the  high  energy  tail  of  the  electron  energy 
distribution,  were  negligible  sources  of  electron  energy  loss. 


A few  predictions  have  also  been  made  for  mixtures  of  h^/Ar. 

The  important  additional  process  in  these  mixtures  is  the  momentum 
transfer  collisions  of  electrons  with  the  inert  gas.  The  cross  section 
data  for  argon  was  taken  from  Frost  and  Phelps.  On  comparing  pure 
gas  calculations  with  mixture  predictions,  it  was  found  that  at  the  same 
E/Nj^2>  the  electron  distribution  of  the  mixture  has  a lower  average 
electron  energy  and  a more  non-Maxwellian  distribution  than  that  of  the 
pure  gas.  In  particular,  the  higher  energy  portion  of  the  electron  energy 
distribution  of  the  mixture  is  decreased  relative  to  that  of  the  pure  gas. 
The  magnitude  of  this  effect  depends  upon  the  ratio  of  diatomic  to  inert 
gas  and  the  value  of  E/N^  . 


This  behavior  is  intimately  connected  with  the  deep  Ramsauer 
minimum  occuring  in  the  low  energy  region  of  the  argon  momentum  trans- 
fer cross  section.  At  values  of  E/Npj^  corresponding  to  low  characteristic 
energies,  i.  e.  , < 1 eV,  the  effect  of  argon  on  the  mixture  transport 
properties  is  minimized  because  of  the  low  cross  section  for  momentum 
transfer  at  these  low  electron  energies.  This  Ar  momentum  transfer  cross 
section  rises  rapidly  with  increasing  electron  energy,  and,  thus,  will  effect 
the  high  energy  tail  of  the  electron  distribution. 


The  magnitude  of  this  effect  is  shown  in  Table  1-1  where  predicted 
discharge  properties  are  shown  for  mixtures  of  50%  H^/bO^o  Ar,  15%  H^/ 
85%  Ar,  and  pure  All  calculations  were  performed  for  a constant 

E/NH?  of  8 x 10~16  v-cm^.  The  most  interesting  comparison  is  between 
the  characteristic  and  average  electron  energies.  In  the  case  of  a 
Maxwellian  distribution  these  quantities  would  be  representative  of  kTe 
and  3/2  kTg  respectively,  where  Te  is  the  electron  temperature^.  As 
can  be  seen,  these  two  quantities  approach  each  other  with  increasing  argon 
concentration,  i.e.  the  electron  energy  distribution  becomes  more  non- 
Maxwellian  at  the  expense  of  the  high  energy  tail  of  the  distribution.  As 
can  be  seen  from  Table  1-1  this  "truncation"  of  the  electron  energy 
distribution  has  a dramatic  effect  on  the  discharge  energy  allocation, 
channeling  energy  away  from  dissociation  into  vibrational  excitation.  It 
is  to  be  noted  however,  that  the  vibrational  excitation  rate  constant  does 
not  vary  significantly  with  mixture  ratio.  This  is  because  the  drift 
velocity  decreases  with  increasing  argon  concentration. 
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TABLE  1-1 

EFFECT  OF  Ar  ON  TRANSPORT  PROPERTIES  OF  H2 
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(all  calculations  performed  for  E/N  =8x10  v-cm  ) 

H 


GAS  MIX 

100%  h2 

50%  H 

2 

50%  Ar 

15%  H2 
85%  Ar 

Characteristic  Energy,  eV 

2.  63 

2.  56 

2.  1 

Average  Energy,  eV 

3.  83 

3.  45 

2.4 

Drift  Velocity,  cm/sec 

8. 5 x 106 

5. 3 x 106 

3. 2 x 106 

Fractional  Discharge  Energy 
Lost  to  Vibrational  Excitation 

0.29 

0.  46 

0.  68 

Fractional  Discharge  Energy 
Lost  to  Dissociation 

0.  55 

0.  32 

0.  02 

Vibrational  Excitation  Rate 
Constant  v = 0 -*  1,  cm^/sec 

-9 

3. 3 x 10 

-9 

3.  3 x 10 

-9 

3x10 
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As  alluded  to  earlier,  this  effect  of  Ar  becomes  less  pronounced  with 
decreasing  E/N^,  because  the  fraction  of  electrons  which  would  fall  into 
the  higher  energy  region  dominated  by  the  argon  momentum  transfer  cross 
section  decreases  with  decreasing  characteristic  electron  energy.  The 
importance  of  the  argon  Ramsauer  minimum  in  the  transport  behavior  of 
H2/Ar  discharges  had  previously  been  discussed  by  Englehardt  and  Phelps.15 
Although  the  H2  cross  sections  employed  in  the  present  analysis  are 
different  than  those  used  in  Ref.  15,  the  phenomenology  is  the  same. 

Predictions  of  the  energy  deposition  and  excitation  rate  constants 
for  pure  D2  are  shown  as  a function  of  E/N  in  Figs.  1-8  and  1-9.  The 
results  are  quite  analogous  to  those  for  H2,  except  that  the  optimum 
discharge  conditions  for  vibrational  excitation  are  at  slightly  lower  E/N 
and  the  maximum  amount  of  energy  transferred  to  vibration  is  slightly 
higher.  No  predictions  have  been  provided  for  D2/Ar  mixtures;  however, 
it  is  expected  that  the  mixture  behavior  would  be  similar  to  that  observed 
for  H2/Ar. 

Discharge  calculations  have  been  performed  for  D2/He  mixtures.* 

A composite  of  the  momentum  transfer  cross  sections  for  He  determined 
by  Frost  and  Phelps14  and  Crompton  et  al,  16  was  employed  in  these 
calculations.  The  predicted  characteristic  energy  and  drift  velocity  for 
mixtures  of  100%  D2,  40%  D2/60%  He,  and  20%  D2/80%  He  are  shown  in 
I- 1 0 as  a function  of  E/Nq  . As  can  be  seen,  both  characteristic 
energy  and  drift  velocity  decrease  monotonically  with  increasing  helium 
concentration.  Unlike  argon,  helium  exhibits  a relatively  smooth  momentum 
transfer  cross  section  with  no  Ramsauer  minimum,  thus  the  electron 
energy  distribution  is  not  truncated.  Indeed,  the  ratio  of  average  to 
characteristic  electron  energy  approaches  1.5,  the  Maxwellian  distribution 
value,  with  increasing  He  concentration. 

The  predicted  discharge  energy  deposition  and  excitation  rate 
constants  for  D2/He  mixtures  of  20,  30  and  40%  are  shown  in  Figs.  I- 1 1 
and  1-12.  These  latter  two  figures  are  plotted  vs  E/Nt  (Nt  = total  gas 
density)  rather  than  E/Nj-j^.  As  can  be  seen,  the  peak  efficiency  for 
vibrational  excitation  is  essentially  the  same  for  all  mixtures.  Thus,  the 
addition  of  helium  would  not  be  kinetically  detrimental  to  laser  performance 
as  far  as  the  discharge  excitation  is  concerned. 


* 

(At  the  time  of  this  analysis  it  was  thought  that  He  would  be  the  carrier 
gas  in  the  D2/HC1  EDL  experiment.  As  is  discussed  in  the  following 
sections  it  has  since  been  found  that  Ar  would  be  more  appropriate  because 
of  fluid  dynamic  mixing  considerations.) 
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Internal  Modes  of  Pure 


VIBRATION 


Comparison  of  Percent  Discharge  Energy  Deposition 
into  Internal  Models  of  D,  for  Various  D./He  Mixtures. 
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Fig.  1-12  Comparison  of  Excitation  Rate  Constants  of  for  Various 
D2/He  Mixtures. 


4.  Plenum  Conditions 


From  the  rate  constants  for  vibrational  excitation  such  as  shown  in 
Figs.  1-4,  1-7,  1-9  and  1-12  the  actual  vibrational  populations  that  result 
from  a given  experimental  discharge  condition  can  be  calculated.  A 
specific  calculation  will  be  carried  out  to  illustrate  the  technique.  The 
chosen  experimental  conditions  are  listed  in  Table  1-2. 

It  will  be  assumed  that  the  discharge  is  uniform  and  fills  the 
plenum.  Then,  the  average  residence  time,  T , for  a molecule  in  the 
plenum  is:  r 


T 

r 


where  u^  is  the  average  gas  velocity.  From  mass  conservation: 


(I-  10) 


m 

UD  p A ’ 


(I-  ID 


Hence,  for  the  given  experimental  conditions,  T =1.8  msec. 

r 

For  an  E/N^  = 0.  8 x 10"^  V-cm^  and  a mixture  of  D2:He  = 0.  33:0.  67, 
the  discharge  calculations  yield  a drift  velocity,  u = 1.9  x 10^  cm/sec. 

The  electron  density,  n^,  may  be  determined  as  follows: 


n 

e 


_I 

u A 
e 


= 5 x 


10  3 

10  electrons/cm 


(1-12) 


Calculated  by  the  code  are  the  rate  constants  for  excitation  into 
the  first  three  vibrational  levels,  and  the  rate  of  transfer  of  energy  into 
rotation  and  translation.  (It  is  assumed  that  the  rotational  and  translational 
modes  are  coupled).  Given  a rate  constant  for  vibrational  excitation,  k., 
the  population  of  the  ith  vibrational  level,  n.,  at  the  end  of  the  discharge 
(plenum)  is 


n.  = k.  n 
l i e 


T 

r 


d-13) 
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TABLE  1-2 


NOMINAL  EXPERIMENTAL  DISCHARGE  CONDITIONS 


Length  of  Discharge  = L = 20  cm 

2 2 

Cross  Sectional  Area  of  Discharge  = A = 1. 9 x 28  cm  - 53.  2 cm 

-17  2 

E/Nt  ~ 8x10  V-cm 


Discharge  Current  = 8 amps 


= Total  Gas  Density  in  particles/cm  of  D^/He  Mixture 


Plenum  Pressure  = 53  torr 


Flow  Rates  of  Gases: 


m = 0.  82  moles/sec 
He 


m =0.41  moles/sec 
D2 


Total  Gas  Density  in  Plenum  (53  torr  at  400  K)  = o = 2.  13  x 10  ^ moles/ 


For  the  particular  experimental  conditions  the  vibrational  populations  that 
result  are: 


,„16  / 3 

n^  =4x10  part/cm 

, , , _ 1 5 . 3 

n^  = 1.3x10  part/cm 

, ,„14  - 3 

n^  = 1 x 10  part/cm 

The  estimate  of  rotational/translational  heating  is  AT  ~ 20°K. 

The  above  calculation  neglects  V -*  T deactivation  from  D^.  We 
estimate  a relaxation  time  of  Ty-p  for  D£  of  ~ 100  msec  for  the  plenum 
conditions.  Clearly,  Tyj  » Tr  and  V -*  T processes  can  be  neglected. 
However,  the  intramode  V -*  V exchange  time  for  D£  is  estimated  to  be 
~ 1 msec  for  the  given  conditions,  and  considerable  rearrangement  of 
quanta  can  occur  within  the  while  the  gas  flows  through  the  plenum. 
Since  the  excitation  rates  to  vibrational  levels  of  higher  than  v = 3 
are  not  known,  and  since  the  V -*  V rates  are  fast,  we  have  assumed  that 
the  distribution  found  in  the  discharge  will  probably  approach  a Boltzmann 
distribution  and  can  be  characterized  by  an  approximate  vibrational 
temperature: 


Ty  = eD  /In  (nQ/n1)  * 1700  K 


(1-14) 


where  9j}  *s  t^ie  characteristlc  vibrational  spacing  for  D2  = 4300  K. 

A more  exact  calculation  would  have  to  include  the  combined  finite 
kinetics  of  electron  excitation  together  with  finite  vibrational  exchange 
rates  as  the  D2/He  mixture  flows  through  the  plenum.  This  type  of 
calculation  is  possible  by  a modification  of  the  present  code,  but  was  not 
pursued  in  the  present  contract. 
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C.  D2/HC1  laser  modeling 

1 . Assessment  of  Vibrational  Energy  Transfer  Rates 

In  order  to  perform  model  calculations  for  the  D2/HCI  laser  it 
is  necessary  to  determine  the  kinetic  rates  for  vibration->trans  lation 
(V-*T)  and  vibration-»vibration  (V-»V)  energy  transfer  processes  for  this 
molecular  system.  Considerable  experimental  data  is  available  for  the 
> HCl,  and  D2/HCI  mixtures,  and  this  data  base  has  been  reviewed  to 
determine  the  best  kinetic  rate  package.  Almost  all  of  these  data  provide 
kinetic  information  on  the  rates  for  the  lowest  vibrational  levels  (v  = 1). 
Therefore,  a combination  of  empiricism  and  theory  has  been  used  to 
extend  the  rates  to  the  upper  vibrational  levels  for  the  anharmonic 
oscillator  model  (A.  H.O.  ). 

The  mechanism  considered  for  vibrational  energy  exchange  in 
the  D2/HCI  system  is  shown  in  Table  1-3.  It  can  be  seen  that  this  is  a 
complex  mechanism  with  many  rates. 

a.  V-»T  Processes 

Considerable  experimental  data  exist  for  the  V-»T  processes  for 
the  D2/D2  and  HC1/HC1  interactions  for  v = 1.  For  D2,  data  exist  at  high 
temperatures  (1500  - 3000°K)  from  shock  tube  experiments  of  Kiefer  and 
Lutz^?  anci  Bird  and  Breshears.  18  At  room  temperature  a number  of 
measurements^"^  have  been  performed.  Quite  recently  Lukas  ik  and 
Ducuing^  have  performed  a Raman  excitation  technique  to  measure  the 
V-»T  relaxation  of  D2  over  the  temperature  range  of  50->400°K.  The 
totality  of  this  data  provide  a consistent  set  of  measurements  and  define 
the  V-*T  rate  constant  and  its  temperature  dependence  to  better  than  a 
factor  of  two.  A curve  fit  to  these  data  is  used  as  the  rate  for  the 
process: 


D2  (v  = 1)  + D2  -»  D2  (v  = 0)  + D2  . (1-15) 

For  the  HC1/HC1  V-»T  process  data  also  exist  over  a wide  range 
of  temperature.^-  ^ Figure  1-13  summarizes  these  data  and  shows  the 
resulting  curve  fit.  It  is  seen  that  this  rate  exhibits  a minimum  around 
room  temperature.  This  phenomenon  is  typical  of  the  hydrogen  halides 
and  is  presumably  due  to  the  dominance  of  attractive  forces  in  the 
molecular  interaction  at  lower  temperatures. 


■I 
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TABLE  1-3 


A NHARMQNIC  OSCILLATOR  KINETIC  MECHANISM 
FOR  D2/HC1 


V -»  T PROCESSES: 

A (v  = 1)  + B 
A (v  = 2)  + B -, 

A (v  = n)  + B -» 
(4  combinations: 

V ->  V PROCESSES: 


A (v  = 0)  + B 
A (v  = 1)  + B 

A (v  = n-1)  + B 
A/A,  A/B,  B/B,  B/A) 


A 

A 

A 

A 

(A 


INTRAMODE 


(V  = 1) 

+ A 

(v  = 1) 

-> 

A(v  = 2)  + A(v  = 0) 

(V  = 2) 

+ A 

(v  = 1) 

-4 

• 

A (v  = 3)  + A (v  = 0) 

(v  = n) 

+ A 

(V  = 1) 

• 

• 

-4 

• 

A (v  - n + 1)  + A (v  = 0) 

(v  = n) 

+ A 

(v  = m) 

• 

• 

-4 

A (v  = n + 1)  + A (v  = m ■ 

= D^  or  HC1,  two 

sets  of  rates) 

1) 


A (v  = 1) 

A (v  = n) 
(A  = D2, 


INTERMODE 


+ B (v  = 0) 

-»  A 

+ B (v  = m) 

L A 

B = HC1) 

(v  = 0)  + B (v  = 1) 

(v  = n-1 ) + B (v  = m + 1) 
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mm 


300 


200 


TEMPERATURE, °K 

o ZITTEL  & MOORE,  1973 
□ HOPKINS,  CHEN,  a SHARMA,  1973 
--BRESHE A RS  a BIRD,  1970 
— FIT  TO  DATA 


k = 4.68  x |0’10  exp(~ 76.83  T'l/3 
+ 3.13  x |0"16  exp (27.91  T'l/3) 


n*j 


(TEMPERATURE, °K) 


V-*T  Data  for  HCl/HCl  Collisions  as  a Function  of 
Temperature  and  Assumed  Curve  Fit. 


There  is  considerably  less  data  for  the  mixed  V -*  T processes 
D2/HCI  and  HCI/D2.  In  fact,  no  measurements  exist  for  D2/HCI. 
Theoretically  HC1  should  be  a less  efficient  collision  partner  than  D2 
in  V -» T deactivation  of  D2.  We  will  assume  that  kD2/HCl=0-2  kD2^D  • 
This  process  is  obviously  not  important  in  the  mechanism  unless  2 

[HC1]  > [D2]  which  is  not  the  optimum  laser  condition.  For  HCI/D2 
there  is  a measurement  by  Hopkins  and  Chen1^  which  indicates  that  D2 
is  about  0.025  as  efficient  as  HC1  itself.  Therefore,  we  have  assumed 
that  kHCL/D2  = 0.  025  kHCl/HCl* 

b.  V-»V  Processes 


For  the  D2/HCI  V ->V  transfer  rate,  experimental  data24  exist  over 
the  temperature  range  of  196  - 342  K.  These  measurements  show  an 
inverse  temperature  dependence  which  has  been  explained  theoretically2^ 
as  due  to  long  range  intermolecular  forces.  Based  on  this  theory  we 
predict  that  k^^y^^cc  T-^  /2.  At  sufficiently  high  temperature  short 


range  forces  may  begin  to  dominate  resulting  in  a different  temperature 
dependence.  However,  for  the  temperatures  of  interest  in  this  program 
(T  < 1000  K)  the  long  range  interaction  should  dominate. 


c.  Anharmonic  Oscillator  Rates 


The  rationale  for  choosing  the  rate  constants  for  v > 1 for  the 
A.  H.  O.  model  is  as  follows:  For  the  V-T  processes,  scaling  to  upper 

vibrational  levels  was  accomplished  by  making  use  of  the  Landau -Telle r 
selection  rules'^  together  with  the  v - 1 -»  0 rates  already  established. 
This  scaling  is  an  approximation,  since  it  was  initially  derived  for 
harmonic  oscillators. 

In  addition  to  this  Landau-Teller  scaling,  it  is  expected  that  there 
will  be  some  additional  change  in  rate  constant  due  to  the  fact  that  the 
exothe rmicity,  AE,  of  the  process  will  depend  upon  the  particular  vibra- 
tional level(s)  involved.  For  V-T  processes  and  V -V  processes  involving 
large  exothe  rmicity,  both  of  which  primarily  involve  short  range  inter- 
actions, the  /\E  dependence  of  the  probability  is  scaled  from  the  expected 
theoretical  dependence.  2? 


k oc  T 1 exp  £-C  (AE)2/3  T‘1/3  J (1-16) 

The  constant  C is  determined  from  the  temperature  dependence  of  the 
v = 1 process. 
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For  V -V  processes  involving  small  fcE  ( ■s  200  cm"  ),  the  energy 
and  temperature  dependence  is  scaled  according  to  long  range  interaction 
theory : 


k <x  T 


-1 


(1-17) 


where  the  constant  C has  been  determined  from  empirical  correlations 
or  from  experimental  data^9-31  where  available. 


By  this  combination  of  experimental  data,  empirical  correlations, 
and  theory,  the  complete  A.H.  O.  kinetic  package  was  derived  and  is 
summarized  in  Table  1-4.  Due  to  their  complexity,  the  A.  H.  O.  kinetic 
rate  constants  are  obviously  an  approximation  at  this  stage,  and  further 
experimental  and  theoretical  investigation  are  required  to  improve  our 
understanding  of  these  processes  for  the  D^/HCl  molecular  system.  The 
rate  constants  listed  in  Table  1-4  are  those  that  were  usea  in  the  modeling 
discussed  in  this  report. 


Toward  the  end  of  this  program,  additional  data  became  available 
that  would  modify  some  of  the  rate  constants  derived  above.  In  particular, 
under  another  program  at  PSI,  Rosen  et  al^  have  recently  measured  the 
high  temperature  (1000  2000°K)  HCI/H2  V ->  T rate.  Combined  with  the 
lower  temperature  data  of  Bott  and  Cohen, 33  these  results  show  that  the 
HCl/H^  V ->  T process  has  a steeper  temperature  dependence  than  that  for 
HCl/HCl.  Thus,  the  assumption  of  kj-jc^/D^  = 0.  025  k 1 / HC 1 *s  Pr°t,ably 
not  correct,  although  the  difference  in  the  vicinity  of  room  temperature, 
where  the  present  modeling  is  concerned,  is  not  expected  to  be  very  signi- 
ficant. In  addition,  Bott  and  Cohen^3  have  also  measured  the  D2/HCI 
V -*V  rate  over  the  temperature  range  of  300  -750°K.  Their  result  shows 
a less  steep  temperature  dependence  than  the  data  of  Hopkins  et  al.^4 
However,  in  the  vicinity  of  300°K  the  results  are  nearly  identical. 

It  has  not  been  possible  to  evaluate  the  quantitative  significance  of 
these  revised  rate  constants  on  the  laser  calculations  presented  in  this 
report.  For  the  cavity  conditions  of  interest  in  this  study  T £ 300°K. 

At  300°K,  there  is  little  difference  in  the  rate  interpretations.  However, 
at  lower  temperatures,  one  would  expect  that  using  the  more  recent  rates 
would  lead  to  less  HCI/D2  V->T  deactivation  and  perhaps  somewhat  lower 
D2/HCI  V-»V  pumping.  These  effects  are  expected  to  be  small  and  of 
minor  importance  in  the  overall  conclusions. 
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TABLE  1-4 


FOR 


D. 


HC1 


LIST  OF  RATES  FOR  D^/HC1  A.H.O.  MODEL 

3 - 1 

(all  k's  in  units  of  cm  -sec  ) 


uu  ' = 31 1 8.  4 cm 

e 

U)  'x  1 = 64.  09  cm 

e e 

uu  ' - 2uu  'x  ' = 2990.  22  cm' 

e e e 

Vibrational  Levels  -*  u s 10 


AE'  = uu  1 - 2u  uu  1 x ' 
e e e 


uu  = 2989.  74  cm 
e 

uu  x = 52.  05  cm  ^ 

e e 

uu  - 2uu  x = 2885.  64 

e e e 

Vibrational  Levels  -»  v s 15 


1) 


2) 


A E = uu  - 2v  uu  x 
e e e 

D2  (U)  + °2  **  °2  (U_1)  + °2 

U,U_1  9 7/3  1/3, 

k (D2'  D2}  = (u)  2-  9°  x 10  7 exp  [-0.  5963  (AE1)  ' T ' ] 

U-l.U  u,  u- 1 

k (D2’  D2)  = k (D2’  °2)  exp  [-1*4388  (AE‘)  T_1] 

Dz  (u)  + HC1  ->  D (u-1)  + HC1 
u,  u-  1 u,  u-  1 


k (D2>  HC1) 
u-  1,  u 

k (D2>  HC1) 


0.  200  k (D2,  D2) 
u,u-l 

k (D2,  HC1)  exp  [-1.4388  (AE1)  T ] 
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TABLE  I_4  (Cont.  ) 


HC1  (v)  + HC1  -*  HCl(v-l)  + HC1 
v,  v-  1 


v’v_i  ( -10  r 2/3  -1/3T 

k (HC1,  HC1)  = (v)  ^4.  676  x 10  exp  [ - 0.  379  (A  E)  T J 

+ 3.  130  x 10'16  exp(27.  91  T_1/3)| 


V-  1 , V V,  V-  1 _ | 

k (HCl,  HCl)  = k (HC1,  HC1)  exp  [-1.4388  (AE)  T'  ] 


HC1  (v)  + D2  -*  HCl(v-l)  + D2 


v,  v 


- 1 


V,  V-  1 

k (HCl,  D2)  = 0.025  k (HCl,  HCl) 

V-  1 , V V,  V-  1 _ J 

k (HCl,  D ) = k (HCl,  D2)  exp  [ - 1.4388  (AE)  T_  ] 


D,  (u)  + HCl(v-l)  -*  D (u-1)  + HCl  (v) 

u>  u~  1 ( -11-1  -3  l 

k (D  , HCl)  = (u)  (v)  < 8.  89  x 10  T exp  [ (- 1.  633  x 1 0 ) AE" ] £ 

2 v - 1 , v 

AE"  =AE'-AE=UJ'-2utu,x,-w  +2vU)x 

e e e e e e 

u- 1 , u u, u- 1 _ j 

k (D  , HCl)  = k (D  . HCl)  exp  [ -1.4388  (AE")  T_  ] 

2 v , v - 1 v- 1 , v 


HCl  (v)  + HCl  (v1  - 1)  -»  HCl  (v-  1)  + HCl  (v') 
|ae|  = (2v'  - 2v)  oj  e xe 


V,  V-  1 


k (HCl,  HCl)  = (v)  (v')  | 6.  769  x 10"11  T' 1/?exp[-4.  638x  10'3  'Jae|] 

v'-l»vl  io  i/o  0/0  -1/3) 

+ 1.941X  10  r exp[-(0.  495|  AEl  + 40.  246)  T 1 J 


V- 1 , V V, V- 1 

k (HCl,  HCl)  = k (HCl,  HCl)  exp  [ - 1.4388  |ae|  t-  ] 
v',  v'-  1 v'-  1,  v' 
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TABLE  I_4  (Cont.  ) 


7)  D2  (u)  + D2  (u'-l)  -*  D2  (u-1)  + D2  (u>) 

I AE'I  = (2u'  - 2u)  oj  'x  1 

1 ' e e 

U,  U-  1 V,  V-  1 

k (D2’  D2)  = 0,200  k<HC1*  HC1) 
u'-l,u'  v'-l,v' 

= (u)(u'){l.  354xl0'nT'1/2exp[-4.638xl0'3|AE'|  ] 

+ 3.  882  x 10' 13  T1/2  exp  [ -(0.495 1 AE'I 2/3 + 40.  246)T_l/3]  } 
u- 1, u u, u-  1 

k (D2,  D2)  = k (D2,  D2)  exp  [ - 1. 4388 1 A E 1 1 T ] 
u',u'-l  u'-l,u' 


3r  ?■■**■* 
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2.  Characteristic  Times 

Before  embarking  on  a detailed  calculation  of  the  laser  cavity,  it  is 
useful  to  examine  the  characteristic  times  for  various  energy  transfer  re- 
actions and  for  the  major  fluid  mechanical  process  of  mixing.  Since  the 
D2/HCI  EDL  modeling  involves  many  parameters  such  as  area  ratio, 
plenum  pressure,  PQ,  and  temperature,  T0,  gas  composition,  cavity 
temperature,  etc.  some  of  which  are  not  independent  quantities,  a 
complete  optimization  is  very  complex.  Instead,  we  will  carry  out  a 
simple  analysis  concentrating  on  establishing  the  approximate  optimum 
Mach  number  at  the  exit  plane  of  the  nozzle  bank,  Me,  for  the  given 
plenum  conditions,  gas  composition,  and  experimental  dimensions. 

For  simplicity  we  will  assume  that  all  energy  transfer  reactions 
occur  after  the  gas  enters  the  cavity.  For  a given  plenum  condition  and 
gas  composition,  the  flow  is  assumed  to  isentropically  expand  through 
the  nozzle.  Hence,  the  thermodynamic  properties  at  the  nozzle  exit  plane 
can  be  parameterized  by  the  nozzle  expansion  ratio,  or  equivalently  Mg. 

We  consider  the  five  fundamental  vibration  transfer  reactions,  i.e., 


* 1 * 

D2  + HC1  -*  D2  + HC1 

* k2 
°2  + °2  - 2 D2 

D2  + HC1  X D2  + HC1 
HC1  + D2  ki  D2  + HC1 


HC1  + HC1  X 2 HC1 


The  asterisk  indicates  that  the  molecule  is  in  a vibrationally  excited  state. 
In  this  approximate  analysis,  we  consider  the  molecules  in  the  harmonic 
oscillator  limit,  hence  intramode  pumping  of  HC1  (processes  1-2)  is 
omitted.  The  characteristic  times,  t's,  associated  with  the  reactions, 
1-18,  can  be  calculated  assuming  constant  temperature  and  density.  For 
example,  the  characteristic  pumping  time  Tj  ~ l/(kj  [HCl]  ) where 
[HCl]  is  the  number  density  of  HCl.  The  characteristic  time  of  a reaction 
is  equivalent  to  the  vibrational  relaxation  time  for  that  process  at  the 
particular  exit  temperature  and  density  that  result  from  the  expansion. 


In  addition  to  the  reaction  times,  there  are  also  the  characteristic 
flow  time,  Tf  = Lc/ue,  where  Lc  is  the  cavity  length  and  ue  is  the  exit 
velocity  and  the  characteristic  mixing  time,  t^.  For  the  experimental 
conditions,  the  mixing  process  will  probably  be  dominated  by  turbulent 
transport.  Sato34  has  shown  that  the  transition  distance  in  a shear  flow 
is  40  to  5 0 times  the  momentum  thickness  at  the  trailing  edge  in  a low 
turbulence  wind-tunnel,  and  the  size  of  the  nozzle  in  the  NRL  device  is 
about  1 mm  while  Lc  ~ 30  cm.  From  a detailed  calculation  of  turbulent 
mixing  in  a shear  layer35,  the  half -width  of  the  shear  layer  is  determined 
to  be  linear  in  x,  i.e.,  6 ~ 0.07  x.  The  density  ratio  between  the  two 
shear  layers  for  this  calculation  is  comparable  to  the  present  experiment. 
Using  the  above  linear  growth  result,  one  can  now  estimate  the  turbulent 
mixing  time,  Tj^j,  which  is  defined  as  the  time  for  the  shear  layer  to  grow 
to  the  width  of  the  nozzle  at  the  exit  plane. 

For  the  conditions  of  PQ  = 60  torr,  TQ  = 300°K,  and  HCI/D2  = 0.  1/ 
0.9,  the  various  characteristic  times  are  shown  in  Fig.  1-14. 

The  V -V  pumping  time  is  seen  to  be  fast  compared  with  all  V-T  deactiva- 
tion times.  Also  shown  are  Tf  corresponding  to  Lc  ~ 30  cm  and  tm. 

As  one  may  observe,  tm  < Tf  which  is  desirable.  However,  the  pumping 
time  approaches  Tf  with  increasing  Me.  For  this  case,  it  appears  that 
the  optimum  condition  should  be  approximately  1 < Me  < 2.  Under  these 
conditions  the  device  would  be  mixing  limited. 


I 


t 


As  shown  in  Fig.  1-15  - 1-17,  similar  maps  have  been  calculated 
for  PQ  = 1 atm,  TQ  = 300°K  and  mixtures  of  10%,  5%  and  2%  HC1  respec- 
tively. Also  shown  on  Fig  J-15  are  lines  for  effective  exit  or  cavity 
temperature,  Te,  and  area  ratio,  AQ/A*. 

A comparison  between  Figs.  1-14  and  1-15  shows  the  following 
features:  For  the  1 atm.  plenum  pressure  (1-15),  all  the  characteristic 
times  for  the  energy  transfer  processes  are  considerably  smaller  at  the 
same  value  of  Me  due  to  the  higher  cavity  pressure.  However,  the  cavity 
flow  times  are  identical.  In  addition,  we  assume  that  the  turbulent  mixing 
time  is  essentially  independent  of  pressure,  3^  and,  hence,  that  Tm  is  the 
same  for  both  conditions.  The  advantage  of  the  higher  plenum  pressure  is 
that  the  pumping  time  may  be  made  much  smaller  than  Tf  and  Tm  for  con- 
siderably larger  exit  Mach  number.  The  larger  Me  leads  to  lower  cavity 
temperature  which  enhances  gain.  Of  course,  a different  discharge  tech- 
nique would  be  necessary  to  provide  volumetric  excitation  of  the  vibrational 
mode  of  D£  at  this  high,  i.e.,  atmospheric,  pressure. 

The  effect  of  varying  HC1  mole  fraction  is  seen  in  Figs.  1-15,  16, 
and  17.  As  the  HC1  concentration  is  decreased,  the  V-»V  pumping  time 
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Fig.  1-16 


Comparison  of  Characteristic  Times  for  a 5%  HC1  - 95% 
, Mixture.  The  Plenum  Pressure  is  Assumed  to  be  1 
mosphere. 


Fig.  1-17 


Comparison  of  Characteristic  Times  for  a 2%  HC1  - 98% 
Mixture.  The  Plenum  Pressure  is  Assumed  to  be  1 
Atmosphere. 


increases,  and  the  dominant  V -»  T deactivation  process  changes  from 
HC1/HC1  to  HCI/D2.  As  mentioned  previously,  this  analysis  ignores  the 
HC1/HC1  V-»V  pumping.  It  is  difficult  to  estimate  a characteristic  time 
for  this  process,  since  it  depends  both  upon  the  concentration  of  vibra- 
tionally  excited  HC1  and  the  number  of  vibrational  levels  being  considered. 
Even  though  the  rate  for  this  process  is  considerably  faster  than  for  the 
D2/HCI  pumping  process  (Table  I-IV),  it  clearly  adds  to  the  overall  time 
to  pump  HC1  to  a sufficiently  high  level  for  maximum  gain.  Thus,  the 
times  required  to  produce  gain  and  lasing  will  be  considerably  larger  than 
shown  in  Figs.  1-14-17. 

The  major  conclusion  to  be  drawn  from  this  analysis  is  that  for 
the  NRL  experimental  device  designed  to  operate  at  Me  > 2,  insufficient 
flow  time  may  be  available  to  provide  for  optimum  pumping.  A higher 
value  of  P0  would  be  desirable,  but  would  require  a major  change  in  the 
excitation  technique. 


3.  Mixing  Analysis 

The  analysis  presented  in  the  preceeding  section  indicated  that 
the  mixing  time,  Tm  « Tf.  Hence,  to  a good  approximation  the  gases 
mix  instantaneously  at  some  exit  properties  at  the  nozzle  exit  plane. 
This  approximation  is  not  reasonable  if  the  pumping  time  (and  other 
energy  transfer  times)  are  smaller  than  Tm.  For  this  latter  condition, 
the  details  of  the  mixing  clearly  influence  the  vibrational  kinetics. 
However,  examination  of  Figs.  1-14-17  indicates  that  conditions  can  be 
found  for  which  Tpump  ~ Tm  « Tf.  Accordingly,  we  have  developed  a 
simple  mixing  model  which  calculates  the  average  properties  at  the 
nozzle  exit  plane  by  assuming  that  the  secondary  injection  mixes  in- 
stantaneously with  the  main  flow. 


A sketch  of  a nozzle  with  secondary  injection  is  shown  in 
Fig.  1-18.  Plenum  gas  (D£  plus  diluent)  enters  at  station  1 with  pressure 
Pi,  density  p and  velocity  uj.  Cold  HC1  is  injected  at  station  2. 

The  gases  are  assumed  to  be  uniformly  mixed  at  station  3,  the  nozzle 
exit  plane.  The  following  analysis  will  provide  the  average  gas  properties 
at  station  3. 


Assuming  one -dimensional  steady  nozzle  flow,  the  governing  equa- 
tion, for  the  Mach  number  variation  between  stations  1 and  2 can  be 
written  as  ^ 


M2) 


? dT 

(1  + yM  ) — !- 


(1-19) 


(1  - M2) 
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lysis. 


where  V is  the  ratio  of  specific  heats,  and  A is  the  nozzle  cross-sectional 
area.  Any  change  in  plenum  temperature  Ts  is  due  to  V-T  energy  transfer. 
With  an  assumed  D2  vibrational  temperature  distribution  one  can  readily 
calculate  the  Mach  number  distribution  and  consequently,  P and  0,  with  the 
aid  of  the  continuity  equation  and  the  equation  of  state. 

Now,  consider  the  flow  through  the  control  surface  of  Fig.  1-18. 

We  shall  write  the  physical  equations  governing  the  flow  in  such  a way  as 
to  obtain  relations  between  the  properties  at  stations  2 and  3: 

t 1 ' 

Continuity:  puA  = u^A^  + ^2  U2  ^2  (1-20) 

P + P 

, r2  ' 

Momentum:  P 2 (A.,  + cos  cp)  + — (A  -A^  -A  ^ cos  cp)  -PA 


and  AQ  is  the  heat  addition  term  which  results  from  the  vibrational  energy 
transfer  between  stations  2 and  3. 


where  quantities  without  subscript  indicate  the  properties  at  station  3, 

X is  the  mole  fraction,  9 is  the  vibrational  spacing,  Ty  is  the  vibrational 
temperature  of  D2  and  fy  is  the  fraction  of  vibrational  quanta  transfered 
from  D2  to  HC1.  The  above  set  of  equations,  Eq.  1-20  to  1-23,  is 
sufficient  for  determining  the  unknowns  P,  T,  p and  u.  Knowing  the 
properties  of  the  secondary  flow,  one  can  rearrange  the  continuity, 
momentum  and  energy  equation  and  after  some  algebra,  obtain 


MW  Pu  A \ 

Bf  ' 17  " '*  + W2  > °2  “2 
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where  w is  the  mass  flow  rate,  F is  the  molar  flow  rate,  and  Cp  is  the 
specific  heat  at  constant  pressure.  With  a known  nozzle  contour,  A/A^, 
one  can  iterate  on  P.  Assuming  P,  one  can  calculate  u from  Eq.  (1-27) 
and  consequently  T from  Eq.  (1-28).  Finally,  from  Eq.  (1-26),  one  can 
calculate  P which  is  then  compared  with  the  assumed  value.  Iteration 
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continues  until  a converged  solution  is  obtained.  The  properties  of  the 
secondary  flow  can  be  determined  from  the  isentropic  relations  if  we 
know  the  stagnation  pressure  and  temperature,  the  gas  composition,  the 
mass  flow  rate  and  the  exit  temperature.  The  results  of  the  above 
instantaneous  mixing  model,  namely  the  thermodynamic  state  of  the 
mixture  at  station  3,  can  be  used  as  inputs  for  the  D2/HCI  EDL  cavity 
calculations. 

4.  Cavity  Model 

To  model  the  D2/HCI  EDL  in  the  cavity  region  a generalized  laser 
kinetics  code  developed  under  another  program38  Was  used.  This  code 
can  integrate  a large  number  of  kinetic  equations  under  a variety  of  fluid 
dynamic  constraints.  For  the  present  calculations  the  cavity  region  was 
assumed  to  be  a constant  area  duct  with  no  radial  gradients  (one  dimen- 
sional). Heat  addition  to  the  gas  from  V-»T  and  V -*V  processes  is 
accounted  for  throughout  the  flow. 

The  model  assumed  11  vibrational  levels  in  D2  and  16  in  HC1. 

This  rather  large  number  of  levels  was  chosen  to  avoid  "barrier  reflection" 
at  the  top  of  the  vibrational  ladder  as  the  intramode  V-W  pumping  takes 
place.  The  results  suggest  that  a smaller  number  of  levels,  particularly 
in  D2,  could  have  been  used  with  some  small  reduction  in  computer  running 
time.  The  vibrational  kinetics  package  and  spectroscopic  constants  used  in 
this  modeling  are  given  in  Table  1-4. 

As  a function  of  flow  distance  (time),  the  code  calculates  the  popu- 
lations of  the  various  vibrational  levels,  small  signal  gain  (or  absorption) 
on  each  rotational  level  of  each  vibrational  level,  and  fluorescence  emis- 
sion per  particle.  In  addition,  a subroutine  can  be  used  to  simulate  power 
extraction  in  the  gain  equals  loss  approximation.  ^9  The  oscillator 
strengths  and  line  width  cross  section  for  HC1  are  those  developed  for  the 
H2/HCI  GDL  modeling.  ^ 

A calculation  to  simulate  a D2/HCI  EDL  experiment  is  performed 
as  follows:  (i)  Given  the  discharge  parameters  and  plenum  gas  conditions, 

the  Boltzmann  code  is  used  to  determine  the  excitation  rate  constants  and 
electron  drift  velocity  as  functions  of  E/N.  (ii)  From  the  results  of  (i) 
a plenum  calculation  is  performed  as  outlined  in  Section  I.  B.  4.  This 
yields  the  vibrational  distribution  of  D2  and  gas  translational  temperature 
at  the  end  of  the  plenum  or  entrance  to  the  nozzle,  (iiil  Given  the  secondary 
gas  injection  conditions  and  nozzle  geometry,  the  instantaneous  mixing 
model  is  used  to  specify  the  properties  of  the  fully  mixed  gas  at  the  exit 
of  the  nozzle.  Vibrational  relaxation  and  further  heat  addition  to  the  gas 
can  be  accounted  for  as  outlined  in  Section  I.C.3.  (iv)  The  output 


conditions  from  (iii)  are  then  used  as  input  to  the  cavity  model.  In 
principle,  the  various  codes  (i)  - (iv)  can  be  combined  so  that  the  entire 
calculation  would  be  performed  in  sequence.  However,  as  is  usually  the 
case  for  complex  modeling  problems,  it  is  much  more  efficient  to  do 
each  calculation  separately  and  provide  some  evaluation  before  proceeding 
onto  the  next  step. 

To  illustrate  the  cavity  model  a sample  calculation  will  be 
presented.  The  discharge  and  plenum  conditions  are  essentially  those 
discussed  in  Sections  I.B.  3 and  4.  As  will  be  remembered,  this 
yielded  a vibrational  temperature  of  D2  ~ 1700°K.  HC1  at  300°K  is 
injected  at  the  nozzle  and  the  input  conditions  for  the  cavity  calculation 
are  shown  in  Table  1-5. 

Results  of  vibrational  population  distributions  of  HC1  at  various 
distances  downstream  of  the  nozzle  exit  are  shown  in  Fig.  1-19.  This 
figure  shows  the  effect  of  anharmonic  pumping  as  quanta  in  the  lower  vibra- 
tional levels  of  HC1  are  pumped  up  the  ladder  by  the  exothermic  V-*V 
transfer  mechanism,  process  (1-2),  resulting  in  inc reased  population  on 
higher  vibrational  levels.  Because  of  the  complexity  of  the  rate  structure, 
it  is  difficult  to  assess  the  exact  mechanisms  involved.  There  are  also 
effects  of  V-»V  transfer  from  D2  to  HC1  and  of  internal  transfer  in  D2. 

At  the  last  station  shown  (x  = 20  cm),  the  translational  temperature  has 
risen  to  205°K,  so  there  has  also  been  some  small  translational  heating 
due  to  V-*  T deactivation  and  non-resonant  V-»V  transfer  processes. 

The  small  signal  gain,  G0,  resulting  from  these  population  distribu- 
tions is  shown  in  Fig.  1-20.  For  each  vibrational  band  only  the  particular 
transition  of  maximum  gain  has  been  plotted.  It  can  be  seen  that  as  the 
V-»V  pumping  of  HC1  occurs,  the  gain  increases  and  shifts  to  higher  vibra- 
tional levels.  The  absolute  magnitude  of  G0  («»  10'^  - 10'^  cm'*)  is 
clearly  too  small  for  a practical  laser.  This  is  to  be  expected  from  the 
low  initial  D2  vibrational  temperature  (Typ^  = 1700  K).  However,  it  is 
possible  to  calculate  the  spontaneous  radiation  (fluorescence)  from  such 
vibrational  distributions  and  compare  these  results  with  experimental 
data. 


^HClJ  = 0.  029  mole  fraction 

[d2]  = 0.333  '•  " 

|He]  = 0.  638  " 

T = 1700°K 
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P = 0.013  atm 
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FRACTIONAL  VIBRATIONAL  POPULATION,  CHC1  w-1 


2.9%  HC1, 33.3%  D2,63.7%  He,  TV[)  = 1700  °K 
Pe  s .013  atm,  Te- 195  °K  ue=l. I km/sec 


VIBRATIONAL  LEVEL,  V hcl 

Vibrational  Population  Distributions  for  HC1  at  Various 
Distances  Downstream  of  the  Nozzle  Exit.  A Straight 
Dine  on  such  a Figure  would  Indicate  a Boltzmann 
D.stribution. 


<*  <u  S 

in  jC  D 
C H XI 
.2  » 

« 2 

in  _ i 

C N c 

J » O 

k 0 .h 

(h  Z .-ti 

« » s 

g«  £ 

.Z, 

2 £ "d 

5 * g 

<u  ° 

U Z 

S 1 2 
.2  * -2 
k O > 
« Q r 

I*  O ™ 

on0 
^ d u 
2 » .0 


d U 
c c 
bO  o g 

•l-<  ■*-<  zj 

[fl  u 3 
- C 6 

d r"3  X 

£ ^ 5 

in  d ^ 


,01  *9  *N1V9  1VN9IS  11VIMS 


D.  COMPARISON  WITH  EXPERIMENTAL  DATA 


1 . Plenum  Conditions 

A series  of  experiments  were  performed  at  NRL  to  attempt  to 
obtain  an  estimate  of  the  D2  vibrational  temperature  produced  by  the 
electric  discharge  in  the  plenum.  Experiments  were  carried  out  in 
which  C02  was  injected  into  the  discharge  excited  D2/He  stream,  and, 
downstream,  lasing  was  observed  on  the  CCL  band  at  10.  6 By  measuring 
the  gain  on  various  v,  J lines  of  the  C02  10.  o\i  band  with  a probe  laser, 
both  the  vibrational  and  rotational  temperatures  of  the  CO2  were  determined. 
Table  1-6  shows  the  results  for  a series  of  runs. 

The  data  of  Table  1-6  distribute  themselves  into  two  classes.  The 
runs  with  low  cavity  pressure,  i.e.,  6-8  torr,  show  a gradually  increasing 
temperature  throughout  the  cavity.  The  runs  at  high  cavity  pressure,  i.e., 
12-14  torr,  show  a rapid  temperature  jump  with  a subsequent  drop  or 
leveling  off  downstream.  This  latter  behavior  is  believed  to  be  associated 
with  some  choking  of  the  flow  at  high  secondary  (CO 2)  gas  injection. 
Nevertheless,  the  data  of  Table  1-6  indicate  CO2  (001)  vibrational 
temperatures  from  1000~2200°K  with  an  average  value  around  1700  K. 

The  results  of  this  experiment  can  be  used  to  infer  an  approximate 
vibrational  temperature  for  the  D2  in  the  cavity.  Two  different  assump- 
tions can  be  made:  (i)  The  D2  vibrational  level  (v  = 1)  is  energetically  in 

close  energy  resonance  with  C02  (01^1).  Thus,  the  D2  (v  - 1 ) mode  can 
be  assumed  to  be  in  rapid  V-»V  equilibrium  with  C02  (01^1).  If  it  is 
further  assumed  that  the  C02  (01^0)  is  rapidly  quenched,  then  for  every 
quanta  of  D2  (v  = 1)  energy  lost  to  C02,  the  C02  (001)  mode  will  gain~70% 
of  that  energy.  Simple  analysis  indicates  that  under  this  situation,  for  a 
measured  C02  (001)  Tv  of  1700°K,  the  Tv  of  D2  must  be  about  2300°K. 

(ii)  In  the  second  case,  we  assume  that  D2  (v  = 1)  and  C02  (001)  are 
directly  coupled  by  a rapid  V-»V  process.  For  the  gas  conditions  of  a high 
D2  vibrational  temperature  compared  to  the  translational  temperature,  the 
Teare -T rainor^®  relation  may  be  used  to  relate  Tv  for  D2  to  that  measured 
for  C02  (001).  This  latter  analysis  leads  to  a value  of  1100°K  for  Tv  of 
D2  for  a measured  value  of  1700°K  for  C02  (001). 

Several  points  can  be  made  concerning  this  experiment  and  the 
resulting  analysis.  First,  the  fluid  dynamic  problems  encountered  in  the 
experiment  suggest  that  the  measurements  and  their  interpretation  should 
not  be  assigned  high  accuracy.  Nevertheless,  the  device  did  demonstrate 
lasing  for  a large  number  of  gases,  and  the  measurement  of  a value  of 

The  temperature  values  at  12.  1 cm  may  be  high  due  to  some  systematic 
e r ro  r . 
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of  1700°K  for  CO^  (001)  is  reasonable  for  a laser  situation.  Secondly,  the 
two  analyses  give  a wide  range  of  possible  D 2 vibrational  temperatures: 

1 1 00-2300°K.  It  is  noted  that  our  discharge  modeling  calculation  gave  an 
intermediate  value,  ~ 1700°K.  Thirdly,  a reasonable  choice  between 
the  two  chosen  D2/CO2  energy  transfer  mechanisms  cannot  be  made 
at  present  on  the  basis  of  existing  experiment  or  theory.  Finally,  a 
more  complex  analysis  using  a finite  kinetic  analysis  and  the  measured 
rate  constants  for  the  D2/CO2  molecular  system  could  be  performed,  but 
does  not  seem  justified  on  the  basis  of  the  uncertainty  of  the  data. 

In  conclusion,  this  experiment  and  analysis  offers  no  detailed 
verification  of  the  plenum  discharge  modeling.  The  results  show  that 
substantial  excitation  of  the  D2  vibrational  mode  must  occur,  and  this  con- 
clusion is  consistent  with  the  results  of  Sec.  I.  B.  Measurements  on  a 
system  more  amenable  to  detailed  analysis  would  be  needed  to  provide  a 
more  accurate  determination  of  D2  vibrational  temperature. 

2.  Fluid  Dynamics 

In  the  early  stages  of  the  experiment  it  was  recognized  by  NRL  that 
for  certain  values  of  secondary  mass  injection,  the  measured  cavity  pressure 
suddenly  increased,  suggesting  a choked  flow  in  the  cavity.  It  was  also 
believed  that  for  values  of  mass  injection  less  than  these  critical  values, 
the  flow  was  reasonably  well  characterized.  It  was  these  latter  (normal) 
conditions  for  which  the  modeling  was  planned.  However,  as  experimental 
data  became  available  it  was  clear  that,  in  general,  the  measured  cavity 
pressure  did  not  closely  correspond  to  the  expected  values  based  on  known 
plenum  pressure  and  nozzle  area  ratio  even  in  the  non-choked  situation. 

A simple  analysis  illustrates  this  point. 

The  detailed  nozzle  dimensions  for  the  NRL  EDL  are  shown  in 
Fig.  1-21.  The  nozzle  throat  height.  A",  is  nominally  1 mm,  so  that  the 
overall  geometric  area  ratio  is  7.  The  case  chosen  for  analysis  used  a 
primary  flow  of  0.33  D2/O.67  He  molar  mixture  through  the  nozzle  config- 
uration shown  in  Fig.  1-21.  For  this  gas  mixture  the  theoretical  flow  rate 
(in  the  absence  of  boundary  layers)  is  given  by 
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at  the  throat,  where  pQ  (atm)  and  TQ  (°K)  are  the  plenum  pressure  and 
temperature,  respectively.  In  the  discharge  region,  where  the  flow  area 
ft*  7 A*,  the  Mach  number  is  — . 08  and  T/Tq  85 . 998,  p/pQ  * . 995.  Thus, 
a static  pressure  measurement  in  the  discharge  region  yields  pQ  within 
0.  5%.  The  selected  plenum  conditions  are  shown  in  Table  I-VH.  For  pQ  = 53 
torr  and  TQ  = 300°K,  the  mass  flow  is  thus  0.632  gm/sec  per  unit  area  of 
primary  throat,  yielding  4.24  gm/sec  total  flow  (A'  =6.7  cm^),  or  — 1.06 
moles/sec,  comparable  with  observed  flow  rates. 
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For  the  case  of  no  secondary  injection,  if  this  D2-He  mixture  were 
expanded  supersonically  to  an  area  ratio  of  7,  the  downstream  pressure 
would  be  ~ ,009  pQ.  Observed  cavity  pressures  were  in  the  range 
0.2  pQ  > p > 0.  1 pQt  (see  Table  1-7).  This  behavior  would  imply  either 
that  the  flow  had  been  shocked  to  a subsonic  exit  velocity  or  that  nozzle  and 
cavity  wall  boundary  layers  were  thick  enough  to  drastically  reduce  the 
effective  area  of  flow  clcs  e to  the  cavity  centerline.  The  observed  cavity 
pressures  for  this  situation  would  correspond  to  an  effective  area  ratio 
A / A*  in  the  range  1.3  - 1.8,  with  corresponding  cavity  temperatures 

0.  57  T > T > 0.  44  T . 

o o 

For  analysis  of  the  case  with  HC1  injection,  the  instantaneous  mixing 
model  was  used  with  a secondary  flow  of  3%  by  volume  of  HC1  (or  2 7%  mass 
ratio)  introduced  at  an  area  ratio  of  2.32,  or  roughly  3 mm  downstream  of 
the  throat.  At  this  point  the  primary  gas  pressure  should  already  be  less 
than  .06  pQ,  and  the  cavity  pressure  for  a uniformly  mixed  supersonic  flow 
should  be  ~.02  p0.  Both  these  pressure  levels  are  lower  than  the  observed 
cavity  pressures  (see  Table  1-7).  However,  diagnostics  employed  with 
CC>2  injected  as  the  secondary  gas  establish  that  some  of  the  cavity  gas  is 
cold  (~200°K  see  Table  1-6),  so  that  it  must  be  assumed  that  the  entire 
cavity  flow  has  not  been  shocked  to  a subsonic  velocity. 

Various  causes  for  this  poorly  characterized  flow  might  be  listed: 

(i)  nozzle  wall  boundary  layers,  (ii)  cavity  wall  (top  and  bottom)  boundary 
layers,  (iii)  inefficient  secondary  injection  with  associated  incomplete 
mixing,  and  (iv)  poor  diffuser /pump  operation.  Estimates  of  most  of 
these  phenomena  were  made  and  none,  by  itself,  appeared  to  be  sufficient 
to  produce  the  observed  flow  effects.  However,  this  analysis  did  suggest 
several  experimental  changes  which  might  be  implemented  to  improve  the 
fluid  dynamics.  For  example,  cavity  boundary  layer  effects  could  be 
minimized  by  diverging  the  top  and  bottom  walls  slightly  to  account  for  the 
displacement  thickness  of  the  boundary  layers.  An  estimate  of  the  behavior 
of  the  HC1  flow  from  the  secondary  injection  ports  indicated  that  the  HC1 
stream  would  not  penetrate  to  the  nozzle  axis.  Penetration  improvement 
could  be  obtained  by  using  fewer  injection  ports  of  smaller  diameter,  or 
by  using  an  inert  carrier  gas  to  permit  higher  values  of  secondary  gas 
injection. 

Toward  the  end  of  the  present  program,  NRL  was  able  to  improve 
, the  flow  in  the  device  substantially.  This  was  accomplished  by  first 

recognizing  that  the  steam  injector  line  was  limiting  the  pressure  in  the 
device  to  about  2 torr.  More  importantly,  the  top  and  bottom  cavity  walls 
were  wedged  at  an  angle  of  1.7°  each  to  compensate  for  boundary  layer 
, growth.  An  example  of  these  improved  flow  conditions  is  shown  in  Table  1-8. 

**  These  data  are  for  CO/He/ C2H2  flows,  but  it  is  clear  that  substantial 

secondary  injection  produces  no  discernable  choking  effects  as  was  shown 
in  the  earlier  data.  Hence,  (in  retrospect)  there  appears  no  evidence  for 
incomplete  mixing  in  the  NRL  EDGDL. 
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Flow  Conditions  with  and  without  Wedged  Cavit 


Because  of  steam  ejector  vent  valve  (upstream  of  the  laser,  near  the  pump),  the  minimum  pressure  in 
the  device  is  about  2 torr. 


3.  Cavity  Fluorescence  Data 

It  was  decided  that  these  aerodynamic  issues  associated  with  the 
experiment  as  it  first  existed,  precluded  a detailed  comparison  of  the 
results  of  the  analysis  with  the  data.  Therefore,  it  was  decided  to  provide 
some  calculations  which  would  encompass  a range  of  possible  experimental 
operation  and  to  provide  a preliminary  design  of  an  improved  device  to  be 
used  for  further  assessment  of  the  viability  of  the  D2 /HC1  EDL. 

A series  of  measurements  were  performed  by  NRL  on  the  infrared 
fluorescence  from  HC1  using  an  optical  system  which  could  be  translated 
parallel  to  the  flow  during  a run.  These  data  consist  of  records  of  infrared 
intensity  in  selected  bandpasses  versus  distance  downstream  of  the  nozzle. 
Examples  of  these  data  are  shown  in  Fig.  1-22  for  an  optical  bandpass  of 
approximately  3.  0 - 4.  Op..  The  top  trace  (a)  is  illustrative  of  an  unchoked 
flow,  while  the  bottom  trace,  (b)  obtained  at  a higher  HC1  injection  rate,  shows 
an  abrupt  leveling  off  of  the  signal  which  is  not  due  to  optical  or  electronic 
saturation.  This  latter  fluorescence  behavior  was  always  associated  with 
the  abrupt  jump  in  cavity  pressure  and  was  believed  due  to  the  presence  of 
shocks  or  choked  flow  in  the  cavity  as  described  previously. 


* 


Although  it  was  decided  not  to  perform  detailed  comparison  of  these 
data  with  the  calculations,  several  qualitative  conclusions  can  be  made. 

(i)  For  the  more  normal  run  Fig.  1-22  (a),  the  fluorescence  intensity  con- 
tinues to  rise  throughout  the  observation  distance  of  about  12  cm.  More 
recent  fluorescent  data  (DARPA-NRL  Laser  Program  Semiannual  Technical 
Report,  1 Jan.  1975-30  June  1975,  NRL  Memorandum  Report  3217,  1976) 
show  that  in  the  non -choked  situation  the  intensities  actually  do  not  show  a 
leveling  off  until  about  20-25  cm  downstream  of  the  nozzle  exit.  This 
observation  can  be  compared  to  the  calculations  of  Fig.  1-20.  The  conditions 
are  not  identical,  but,  nevertheless,  for  these  low  cavity  pressure,  ~8  torr, 
the  HC1  is  calculated  to  be  at  a low  level  of  excitation  with  small  values  of 
gain  at  distances  of  10  cm  downstream  of  the  nozzle.  Considerably  longer 
flow  lengths  and/or  higher  cavity  pressures  would  be  required  to  provide 
optimum  pumping  and  gain,  (ii)  No  lasing  of  HC1  could  be  obtained.  This 
result  is  also  consistent  with  the  results  of  Fig.  1-20.  (iii)  Observations  at 
longer  wavelength  4.0  - 4.5  ^ showed  much  less  radiation  intensity  than  in 
the  3-4p  bandpass.  Again  this  result  is  consistent  with  the  low  level  of  HC1 
vibrational  excitation  expected.  As  the  gas  flows  down  the  cavity,  HC1/HC1 
V-»V  pumping  should  populate  the  higher  vibrational  levels  of  HC1  and  the 
fluorescence  should  shift  to  longer  wavelength,  (iv)  The  typical  E/N  for  the 
runs  of  Fig.  1-22  is  1 x 10  ^ V-cm  . This  is  close  to  the  optimum  discharge 
condition  for  this  mixture  as  seen  from  Fig.  I - 1 1 . 

These  qualitative  comparisons  of  the  fluorescence  data  combined  with 
the  plenum  comparison  (Sec.  I.  D.  1)  do  indicate  a general  consistency  with 
the  modeling  results.  Excitation  of  the  D2  vibrational  mode  in  the  discharge 
appears  to  occur,  and  it  may  be  possible  to  channel  as  much  as  70%  of  the 
discharge  energy  into  the  vibrational  mode  by  operating  at  optimum  E/N. 

V-* V transfer  from  D2  to  HC1  does  occur  as  evidenced  by  the  HC1  fluorescence. 
However,  to  demonstrate  lasing,  higher  density  operation  is  required. 


-57- 


. *tr. 


l 


E.  DESIGN  OF  NEW  EXPERIMENT 

1.  Concept 

Based  on  the  discussions  of  the  proceeding  section,  it  is  clear  that 
a new  facility  is  necessary  for  a meaningful  experimental  demonstration 
of  the  D2/HCI  EDL  concept.  First,  the  flow  in  the  present  device  cannot 
be  well  characterized.  In  addition,  even  if  the  flow  could  be  improved, 
the  calculations  indicate  that  a higher  density  operation  is  desireable. 

This  is  not  feasible  in  the  existing  apparatus. 

In  the  course  of  the  present  program,  NRL  had  planned  to  replace 
the  pin-type,  ballasted  glow  discharge  with  an  e-beam  sustained  discharge. 
This  latter  technique  is  a recognized  scalable  excitation  method  which  can 
be  used  at  high,  e.  g.  , one  atmosphere,  pressure.  Therefore,  it  was  decided 
to  take  as  a starting  point  for  a new  device  concept  an  atmospheric  pressure 
plenum  excited  by  an  e-beam  sustained  discharge.  Using  the  comparison  of 
the  various  characteristic  times  for  the  major  kinetic  processes  with  the 
fluid  dynamic  constraints  for  the  high  pressure  operation,  the  overall  cavity 
Mach  number  and  nozzle  expansion  ratio  can  be  estimated. 

There  were  insufficient  time  and  resources  available  at  this  point 
in  this  program  to  provide  extensive  discharge  modeling  for  the  e-beam 
sustainer  discharge.  However,  there  is  no  reason  to  believe  that  the 
overall  conclusions  of  the  glow  discharge  modeling,  i.  e.  , that  a large 
fraction  of  the  electron  energy  can  be  channeled  into  the  D2  vibrational 
mode,  will  change  for  an  e-beam  sustained  discharge.  Therefore,  a 
series  of  kinetic  calculations  were  performed  for  various  assumed  values 
of  Tv  of  D2  within  the  range  1000°  ^ Tv  £ 2000°K. 

Finally,  it  has  been  shown  for  the  N2/CO2  mixing  GDL4*  that 
injection  of  the  secondary  gas  at  the  throat  of  the  nozzle  provides  more 
efficient  mixing.  Therefore,  it  was  decided  to  design  a nozzle  for  this 
new  experiment  that  would  provide  for  throat  injection  as  well  as  meet 
the  requirements  for  providing  optimum  cavity  temperature  and  pressures 
as  determined  from  the  kinetics.  A preliminary  design  of  such  a nozzle 
array  was  undertaken. 

The  following  sections  present  the  results  of  this  design  study  and 
can  be  used  to  provide  guidance  in  planning  an  upgraded  experimental 
facility  in  which  to  demonstrate  the  D2/HCI  EDL. 

2.  Kinetic  Modeling 

As  was  discussed  in  Sec.  I.  C.2,  a comparison  of  the  characteristic 
times  for  the  various  vibrational  transfer  processes  with  the  approximate 
mixing  time  and  cavity  flow  time  provides  a method  for  quickly  optimizing 
the  overall  fluid  dynamics  of  the  device.  Various  maps  for  one  atmosphere 


plenum  pressure  and  several  HC1  concentration  were  shown  in  Figs.  1-15  to 
1-17.  These  maps  did  not  consider  any  diluent  in  the  plenum  gas.  An 
additional  map  for  Ar  diluent  is  given  in  Fig.  1-23.  Although  the  bulk  of 
the  measurements  with  the  existing  device  were  done  with  He  diluent,  Ar 
would  be  a better  choice  for  the  new  experiment.  Ar  would  provide  more 
stopping  power  than  He  for  the  high  energy  electrons  of  the  e-beam  in 
the  plenum,  and,  also,  the  Ar/D2  mixture  would  provide  a better  molecular 
weight  match  with  the  HC1  secondary  flow  and  thereby  reduce  dissipative 
effects  in  the  mixing. 

Figure  1-23  shows  that  for  > 2 the  mixing  time  is  less  than  the 
D2/HCI  V-*  V transfer  time  and  both  time  are'  much  shorter  than  the  flow 
time  for  a 30  cm  long  cavity.  For  purposes  of  a design  point  we  have 
chosen  to  select  Me  = 2.  5 to  produce  the  required  cavity  conditions.  This 
requires  a nozzle  of  area  ratio  ~ 3.  3. 

A series  of  full  vibrational  kinetic  calculations  were  performed 
under  the  following  assumptions:  (i)  The  plenum  gas  is  a Ar  mixture 
of  molar  ratio  1/2,  PQ  = 1 atm.,  and  TQ  = 300°K.  Various  D2  vibrational 
temperatures  were  assumed  covering  the  range  of  1000  - 2000°K.  (ii)  For 
each  value  of  D2  Tv,  three  different  amounts  of  HC1  are  injected  at  the 
throat,  and  the  instantaneous  mixing  model  (Sec.  I.  C.  3)  is  used  to  specify 
the  fully  mixed  conditions  at  the  exit  of  the  nozzle.  The  expansion  area 
ratio  is  chosen  to  be  3.  3 to  provide  a nominal  cavity  Mach  number  of  2.  5. 

It  is  further  assumed  that  no  vibrational  kinetics  occur  in  this  mixing 
process,  (iii)  The  fully  mixed  gas  enters  the  cavity,  the  vibrational  kinetics 
occur,  and  small  signal  gains  for  HC1  v,  J transition  are  calculated  as  a 
function  of  flow  distance  through  the  cavity. 

An  example  of  such  a calculation  is  shown  in  Fig.  1-24  for  an  assumed 
D->  T of  1500°K  and  a final  HC1  mole  fraction  of  0.  006  (comparable  to  the 

“ V 

conditions  of  Fig.  1-23).  It  should  be  noted  that  to  preserve  clarity  in  the 
figure,  not  all  the  vibrational  transitions  showing  gain  have  been  plotted. 

Several  points  can  be  made  from  Fig.  1-24.  First,  the  gain  peaks  a few 
cm  downstream  of  the  throat.  This  is  to  be  compared  with  the  results  of 
Fig.  1-20  and  is  due  primarily  to  the  much  higher  cavity  pressure  (45  torr) 
compared  to  the  previous  condition  (~  1.0  torr).  Second,  the  predicted 
small  signal  gains  are  high,  > 1%  cm'*,  and  lasing  should  be  possible  with 
reasonable  values  of  optical  cavity  parameters. 

To  illustrate  the  effect  of  D2  vibrational  temperature,  a typical 
transition  was  selected  (v  = 3 -*  2),  and  the  peak  value  of  gain  for  that 

transition  obtained  in  each  calculation  plotted  against  the  assumed  D2  Tv  j 

for  that  case;  the  result  is  shown  in  Fig.  1-25.  It  can  be  seen  that  even  for 
rather  modest  values  of  Tv  = 1250°K,  reasonable  gains  are  obtained.  From 
the  results  of  the  discharge  modeling  as  well  as  the  D2/CO2  gain  measure- 
ments presented  previously,  we  have  every  reason  to  believe  that  values 
of  D2  Tv  from  1000  - 2000°K  are  obtainable  via  practical  discharge 
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Fig.  1-24  Small  Signal  Gain  for  Selected  Vibrational  Transitions 
of  HC1  as  a Function  of  Distance  Downstream  of  the 
Nozzle. 
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techniques.  Therefore,  the  results  of  this  kinetic  analysis  indicate  that 
a D2/HCI  EDL  operated  at  high  density  using  e-beam  sustainer  discharge 
excitation  would  be  feasible. 

3.  Nozzle  Design 

One  of  the  possible  configurations  for  a D^/HCl  mixing  EDL  would 
be  to  expand  the  D2  through  a two-dimensional  supersonic  nozzle,  while 
the  HC1  is  injected  at  the  throat.  Taran  et  al4^  have  shown  theoretically 
and  experimentally  that  this  configuration  gives  excellent  results  for  the 
N2/CO2  EDL.  It  is  shown  that  satisfactory  mixing  can  be  obtained  with 
an  increase  of  the  freezing  efficiency  over  simple  expansion  nozzles  or 
premixed  flow. 

We  have  chosen  to  design  a nozzle  for  the  next  D2/HCI  EDL  device 
incorporating  throat  injection.  A complete  turbulent  mixing  calculation 
is  beyond  the  scope  of  the  present  program,  but  we  will  use  an  existing 
PSI  laminar  mixing  code  with  appropriate  modifications  as  a design  tool 
for  this  nozzle.  The  turbulent  transport  will  be  modeled  by  the  eddy 
viscosity  approximation. 

There  are  various  successful  eddy  viscosity  models,  each  developed 
for  a different  flow  situation.  ^2,  43  por  pianar>  free  turbulent  mixing 
in  a coflowing  stream,  we  choose  Schetz's44  model  in  order  to  account  for 
the  variable-density  situation  through  the  appropriate  definition  for  the 
mass  flow  defect  (or  excess): 


e = 0.  0036  o 


P u 
e e 


(1-30) 


where  € is  the  eddy  viscosity,  p is  the  density,  u is  the  velocity,  and  subscript 
e indicates  the  external  condition. 

Before  we  apply  this  model  to  our  rase,  it  was  vei  fied  against  two 
turbulent  wake-flow  measurements,  one  oi  which  is  subsonic  and  the  other  super- 
sonic.  Turbulence  measurements  in  the  wake  of  a thin  flat  plate  were 
performed  by  Chevray  and  Kovasznay4^  in  a low-speed  wind  tunnel.  This 
experiment  offers  a critical  test  of  the  present  model,  since  it  corresponds 
to  the  development  of  a turbulent  boundary  layer  with  zero  pressure  gradient 
when  the  wall  shear  is  suddenly  removed  at  the  trailing  edge.  The  results  of 
the  model  comparison  with  the  data  are  shown  in  Fig.  1-26.  The  decay  of 
the  center  line  velocity  defect 
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(1-31) 


Fig.  1-26  Comparison  of  Predictions  of  Schetz's  Model  with 
Subsonic  Data. 
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is  plotted  as  a function  of  x/0  where  0 is  the  momentum  thickness  at  the 
trailing  edge.  The  adequacy  of  Schetz's  eddy  viscosity  model  for  this  low 
speed  case  is  clearly  demonstrated  by  these  results. 

Although  the  two-dimensional  turbulent  wake  has  been  extensively 
studied  at  low  speed,  information  for  supersonic  turbulent  wakes  is  scarce. 
Measurements  with  cylinders  placed  normal  to  a Mach  6 air  stream  have 
been  made  by  McCarthy  and  Kubota4^  and  Behrens.  47  In  these  two  experi- 
ments, the  combination  of  high  Mach  number  and  severe  body  bluntness 
produced  large  gradients  in  the  flow  immediately  adjacent  to  the  wake. 

Under  such  circumstances  the  turbulent  wake  mean  flow  becomes  indistin- 
guishable from  that  of  the  inviscid  wake.  To  restore  the  prominence  of  the 
turbulence  core,  the  inviscid  gradients  must  be  reduced.  This  can  be 
accomplished  by  using  a very  slender  body.  Demetriade s48  has  made 
some  mean-flow  measurements  in  the  supersonic  wake  of  a very  slender, 
two-dimensional  body  in  a continuous  open-circuit  wind  tunnel  at  Mach  3 
with  a stagnation  pressure  of  730  mm  Hg  absolute  and  a stagnation  tempera- 
ture of  38°C.  Under  these  conditions,  the  role  assumed  at  low  speeds  by 
the  cylinder  diameter  as  a gaging  length  of  the  wake  flow  is  assumed  by 
the  wake  drag  thickness.  If  the  body  drag  coefficient,  CD,  is  associated 
with  a lateral  dimension  h,  then  the  wake  drag  per  unit  span  is: 


— Du  h = / p u (u^  - u)  dy  - constant  . (1-32) 


The  constant  quantity  CDh’  called  "the  virtual  body  thickness",  can  be 
obtained  from  this  expression  if  we  assume  that  the  lateral  velocity  distri- 
bution is  similar.  From  the  experimental  data,  the  constant  Cj-,h  is  found 
to  be  0.  00909  cm,  and  this  value  is  used  to  nondimensiortalize  the  axial 
distance.  As  shown  in  Fig.  1-27,  the  rate  of  decay  of  the  centerline 
velocity  defect  is  accurately  predicted  by  Schetz's  eddy  viscosity  model 
for  the  turbulent  transport  in  this  supersonic  flow  situation. 

Having  demonstrated  by  these  comparisons  with  turbulent  flow  data 
that  Schetz's  model  is  reasonably  accurate,  we  can  now  apply  it  to  determine 
the  appropriate  mixing  nozzle  for  the  D2/HC1  EDL.  A sketch  of  the  assumed 
nozzle  configuration  is  shown  in  Fig.  1-28.  The  HC1  injection  port  is  located 
at  the  sonic  point  of  the  nozzle  and  provides  co-axial  injection  of  HC1  or 
HCl/Ar  mixtures.  In  practice,  the  position  of  the  HC1  injection  relative  to 
the  throat  could  be  varied  to  provide  for  some  empirical  optimization  of  the 
mixing.  The  nozzle  throat  height  is  assumed  to  be  1.  0 mm  with  a nominal 
A/A'  of  3.3  to  provide  for  the  required  expansion  to  Me  = 2.5. 

For  this  calculation,  the  laminar  coefficient  of  viscosity,  which  can  be 
readily  calculated  from  the  molecular  properties4*^,  was  determined,  and  was 
found  to  be  small  compared  with  the  eddy  viscosity.  Thus,  our  assumption 
that  shows  that  turbulent  mixing  occurs  appears  to  be  valid. 

-67- 


— — gjftr',...  Bfi __  ■ 


Fig.  1-28  Schematic  of  Mixing  Nozzle  with  Axial  Injection  at  Throat. 

This  Configuration  was  used  in  the  Mixing  Calculation 
Described  in  the  Text. 


We  have  modeled  a case  in  which  the  plenum  flow  contains  approxi- 
mately equal  mole  fractions  of  and  Ar  at  one  atmosphere  pressure  and 
300°K.  The  secondary  flow  consists  of  HCl/Ar  0.  18/0.82  at  1.02  atmospheres 
and  300°K.  The  slight  over  pressure  of  secondary  flow  is  required  to 
maintain  a uniform  static  pressure  at  the  throat  where  both  flows  merge. 

The  calculation  starts  at  the  throat  where  the  Mach  number  is  unity.  The 
velocity  and  density  profiles  at  the  nozzle  exit  plane,  x = 1.  5 cm,  are  given 
in  Figs.  1-29  and  1-30  respectively,  as  a function  of  the  normalized  radial 
distance. 


r = r/rg  , (1-33) 

where  the  nozzle  exit  radius,  rg  = 13.  5 mm.  The  concentration  profiles  for 
and  Ar  are  shown  in  Fig.  1—31  for  various  axial  distances.  It  can  be 
seen  that  the  gases  are  reasonably  well  mixed  at  x ^ 10  cm.  The  mixing 
of  the  HC1  throughout  the  flow  is  shown  in  Fig.  1-32,  where  the  profiles  for 
HC1,  expressed  as  % of  HC1  relative  to  only,  are  given  for  several  down- 
stream stations.  Again,  it  is  seen  that  HC1  is  reasonably  well  mixed  for 
x ^ 10  cm. 

The  conclusion  of  this  analysis  is  that  for  centerline  injection  of 
HC1  at  the  throat  of  the  nozzle,  reasonable  mixing  will  occur  in  a distance 
of  about  10  cm  for  gas  conditions  and  mixture  compositions  which  kinetically 
should  produce  high  gains  and  lasing  at  comparable  distances. 


Fig.  1-29  Velocity  Profile  at  Exit  of  Nozzle  as  Determined  from 
Mixing  Calculation. 
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Density  Profile  at  Exit  of  Nozzle  as  Determined  from 
Mixing  Calculation. 


F.  CONCLUSION  AND  SUMMARY 


An  initial  goal  of  this  program  was  to  develope  an  analytic  capabi- 

o'f  threm°d  the  Dz/HC1  EDL  exPeriment«  This  code  development  consisted 
of  three  major  separate  tasks:  (i)  a model  of  the  discharge  or  plenum 

conditions,  (n)  a technique  to  handle  secondary  injection  and  subsequent 
ga  mixing,  and  (in)  a vibrational  kinetic  code  to  calculate  gain  onvibra- 
tional/ rotation  transitions  in  the  cavity.  These  codes  were  developed  and 
successfully  operated  under  this  program. 


To  specify  the  discharge  conditions,  an  existing  Boltzmann  trans- 
port code  was  modified  and  used.  State-of-the-art  cross  sections  for  electron 
impact  excitation  of  H2  were  determined  by  a careful  review  of  available 
data  and  theory,  and  by  comparison  of  predicted  transport  properties  of 
H2  discharges  with  measurement.  A cross  section  data  package  for  D 
was  determined  primarily  by  analogy  to  the  data  for  H2.  With  these  2 
data  packages,  the  capability  exists  to  calculate  the  transport  properties, 
excitation  rate  constants,  and  energy  deposition  into  the  various  internal' 
modes  for  a wide  range  of  experimental  discharge  conditions  in  pure  H, 
and  D2  and  in  mixtures  of  these  gases  with  He  or  Ar.  2 


The  results  of  a discharge  calculation  for  conditions  comparable 
to  those  used  in  the  NRL  EDL  device  indicated  that  a D?  vibrational 
temperature  of  about  1700°K  was  expected.  No  direct  experimental 
determination  of  this  temperature  was  possible.  However,  approximate 
D2  vibrational  temperatures  could  be  inferred  from  D2/COz  laser  gain 
experiments  and  indicated  a range  of  temperatures  from  1100  - 2300°K 
It  is  undoubtedly  a coincidence  that  the  median  of  this  temperature  range 
agrees  with  the  prediction. 


An  approximate  mixing  analysis  was  developed  to  provide  for 
injection  of  HC1  within  the  nozzle  and  the  instantaneous  mixing  of  the  gas 
with  the  calculation  of  new,  appropriate  flow  conditions  at  the  nozzle  exit’ 
with  the  specification  of  any  heat  input  to  the  gas  due  to  vibrational  kinetics. 
This  model  seemed  appropriate  for  the  NRL  EDL,  since  it  was  estimated 
that  the  mixing  time  was  probably  short  compared  to  the  characteristic 
time  for  V-*V  pumping  of  HC1  by  D2.  This  mixing  code  uses  as  input  the 
plenum  conditions  (as  specified  by  the  Boltzmann  code)  and  calculates  the 
input  conditions  for  the  cavity  analysis. 

A complete  set  of  vibrational  energy  transfer  rate  constants  for  the 
D2/HC1  molecular  system  was  developed  from  a combination  of  experi- 
mental data,  theory,  and  empirical  correlation.  These  rate  constants 
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were  used  with  an  existing  anharmonic  oscillator  laser  code  to  calculate 
HC1  vibrational  populations  and  small  signal  gain  as  the  gas  flowed  down 
the  cavity.  In  addition,  infrared  fluorescence  per  molecule  and  simulated 
power  extraction  are  other  available  options  from  such  a calculation. 

The  second  major  goal  of  this  research  effort  was  to  use  the 
predicted  results  of  the  analytic  model  to  compare  with  experimental  data 
and  to  guide  the  experimental  effort.  An  early  calculation  performed  for 
appropriate  device  conditions  indicated  that  at  the  low  gas  densities  expected 
in  the  cavity,  the  gain  was  not  expected  to  peak  until  well  beyond  the  avail- 
able, accessible  viewing  area.  This  prediction  was  qualitatively  consistent 
with  the  data;  infrared  fluorescence  (observed  primarily  from  the  lower 
levels  of  HC1)  increased  monotonically  throughout  the  observation  distance 
of  about  2 3 cm  and  no  lasing  could  be  obtained. 

More  quantitative  comparisons  of  the  results  of  the  calculations 
with  data  were  inhibited  by  the  realization  that  the  fluid  dynamic  flow  in 
the  NR  L EDL  device  as  it  then  existed  was  difficult  to  characterize.  It 
was  determined  that  the  cavity  static  pressures  were  much  higher  than 
expected  on  the  basis  of  simple  inviscid  calculations.  The  reason  for  this 
discrepancy  was  eventually  determined  to  be  a combination  of  steam  ejector 
and  cavity  wall  boundary  layer  effects. 

Because  of  these  flow  problems,  a decision  was  made  late  in  the 
program  to  use  the  remaining  resources  to  conceptually  design  a better 
experiment.  As  a starting  point,  it  was  decided  to  consider  an  e-beam 
sustainer  discharge  operating  at  one  atmosphere  plenum  pressure  as  the 
laser  excitation  source.  By  operating  at  a higher  plenum  pressure,  higher 
cavity  density  is  attainable.  Although  detailed  discharge  calculations  were 
not  performed  for  these  higher  pressures,  it  is  expected  that  D2  vibrational 
temperatures  comparable  to  the  glow  discharge  results  should  be  attainable. 
A series  of  cavity  kinetic  calculations  were  performed  for  a range  of  D£ 

T 's  from  1000  - 2000°K.  These  results  indicate  the  high  gains,  i.  e. , 

>1%  cm-1,  can  be  achieved  a few  cm  downstream  of  the  nozzle. 

In  addition,  a turbulent  mixing  calculation  was  performed  to  assess 
axial  injection  of  secondary  gas  at  the  nozzle  throat.  It  was  determined 
that  the  gases  should  be  reasonably  well  mixed  10  cm  downstream  of  the 
injection  position.  Obviously,  further  iteration  can  be  performed  to 
provide  the  optimum  device  conditions,  where  the  mixing  length  and  HC1 
pumping  length  are  comparable. 


t. 
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The  conclusions  of  this  conceptual  design  study  are  that  the  higher 
density  operation  is  desirable  and  that  a demonstration  of  a D2/HCI  EDL 
should  be  possible.  The  preliminary  operating  conditions,  nozzle  expan- 
sion ratio,  and  injection  technique  are  determined.  This  information  can 
be  used  to  develop  a detailed  experimental  design. 
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II.  ANALYTICAL  MODELING  OF  A PULSED  ELECTRIC 
DISCHARGE  HC1  LASER 


A.  INTRODUCTION 

This  section  of  this  report  provides  a description  of  research 
performed  at  PSI  in  support  of  the  NRL  pulsed  electrical  discharge  HCl 
infrared  laser  experiment.  This  latter  program  involved  studies  of  the 
electrical  excitation  of  mixtures  of  HCl/Ar  with  the  goal  of  determining  the 
feasibility  of  developing  an  efficient  laser  operating  on  HCl  vibration/ 
rotation  transitions.  In  this  study  the  electrical  excitation  was  provided 
by  an  electron  beam  sustained  discharge  utilizing  the  NRL  Maxwell  cold 
cathode  e-beam  sustainer  device  which  has  the  capacity  of  operating  at 
voltages  of  300  kv  with  beam  current  densities  of  ~ 10  amps /cm  . The 
laser  test  chamber  in  these  experiments  had  an  active  length  of  100  cm, 
height  and  depth  of  10  cm,  and  a mirror/coupling  configuration  designed 
to  allow  observation  of  lasing  at  small  signal  gains  of  order  10  ^ cm 

The  primary  goal  of  the  PSI  research  effort  was  to  develop  an 
electric  discharge  kin  -ic  model  for  HCl  to  be  used  to  provide  predictions 
of  system  gain  and  optimum  operating  conditions.  A detailed  discussion  of 
this  electron  kinetic  model  is  presented  in  Sec.  Il-B  and  specific  system 
gain  predictions  and  comparison  with  dat"’  are  discussed  in  Sec.  II-C. 
Lastly,  an  analysis  of  infrared  fluorescence  measurements,  performed 
early  in  the  experimental  program,  is  given  in  Sec.  Il-D. 
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ELECTRIC  DISCHARGE  KINETICS 


It  is  well  known  that  when  a partially  ionized  gas  is  subjected  to  an 
electric  field  the  resulting  steady  state  electron  velocity  distribution  need 
not  be  Maxwellian.  Indeed  the  resulting  velocity  (energy)  distribution  is 
determined  by  the  competition  between  the  accelerating  electric  field  and 
elastic  and  inelastic  collisions  between  the  electrons  and  bath  gas  molecules. 
The  characteristics  of  such  discharges  are  generally  defined  in  terms  of  the 
bulk  transport  properties  for  given  values  of  the  Townsend  parameter  E /N, 
where  E is  the  applied  voltage  per  unit  length  and  N is  the  gas  number 
density.  The  two  most  commonly  used  transport  properties  are  the  electron 
drift  velocity  Vq  and  the  characteristic  electron  energy  . In  the  case  of 
a Maxwellian  velocity  this  latter  quantity  is  the  equivalent  of  k Te  where  k 
is  Boltzmann's  constant  and  Tg  is  the  electron  temperature. 

In  order  to  specify  the  excitation  rate  constants  for  specific  electron/ 
molecule  interactions  occuring  in  the  discharge,  which  is  the  goal  of  the 
present  analysis,  the  detailed  electron  velocity  distribution  is  required. 
Although  this  distribution  cannot  be  readily  deduced  expe  rimentally,  it  can 
be  predicted  theoretically  by  solving  the  Boltzmann  transport  equation. 

This  analytical  technique,  however,  requires  specification  of  the  electron 
energy  dependent  cross  sections  for  all  of  the  important  elastic  and  in- 
elastic electron/molecule  processes  which  occur  in  the  discharge. 
Unfortunately,  such  detailed  data  is  not  commonly  available  for  most  gases. 

If  the  transport  properties  of  a gas  are  known  as  a function  of  the 
Townsend  parameter,  then  unknown  cross  sections  may  be  deduced  by  an 
inverse  process,  e.g.,  ref.  50.  In  this  method,  the  Boltzmann  transport 
equation  is  solved  with  a trial  set  of  cross  sections  in  order  to  provide 
predicted  bulk  transport  properties  as  a function  of  the  Townsend  parameter. 
The  predictions  are  then  compared  to  the  data  and  the  trial  cross  sections 
are  adjusted  in  a systematic  manner  in  order  to  improve  the  agreement 
between  prediction  and  data.  This  process  is  repeated  iteratively  until  the 
difference  between  data  and  predictions  is  within  the  uncertainty  of  the  data. 
This  technique  cannot  provide  a unique  set  of  cross  sections,  and,  therefore, 
considerable  insight  is  required  in  choosing  the  trial  cross  sections.  The 
selection  process  can  be  simplified  somewhat  if  the  analysis  is  limited  to  a 
range  of  Townsend  parameter  where  only  one  or  two  inelastic  processes  are 
important. 


In  the  case  of  HC1,  the  data  base  is  relatively  sparse.  Only  a few  of 
the  inelastic  processes  of  interest  in  HC1  have  been  studied,  and,  furthermore, 
the  only  measurement  of  the  bulk  transport  properties  of  HC1  discharges  is 
the  early  work  of  Bailey  and  Duncanson.  5 1 ’ 52  This  latter  study  was  performed 
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in  pure  HC1  over  the  E/N  range  of  3 - 12  x 10'^  V-cm^,  and  the  observed 
values  of  drift  velocity  and  characteristic  energy  are  listed  in  Table  11-1. 

In  the  present  analysis  these  data  will  be  used  in  conjunction  with  the  avail- 
able cross  section  information  to  deduce  a self-consistent  set  of  cross 
sections  for  HC1. 

There  is  one  caveat  which  should  be  mentioned  before  proceeding 
with  this  analysis.  Although  Bailey  and  Duncans  on1  s ^ ^ work  is  qualitatively 
valuable,  their  measured  transport  properties  may  not  be  as  accurate  as  one 
might  desire.  In  particular,  the  drift  velocities  were  determined  through 
use  of  a magnetic  field,  and  the  difference  between  electric  and  magnetic 
drift  velocity  was  not  appreciated  at  that  time.  Indeed,  the  actual  drift 
velocity  can  be  determined  from  their  raw  data  only  if  a number  of  restrictive 

C Q 

assumptions  are  made.  In  the  one  instance  where  modern  measurements 
can  be  compared  with  measurements  performed  by  Baileys'  group, the 
agreement  is  poor,  ^ suggesting  that  the  error  associated  with  Baileys' 
experimental  technique  is  large.  Thus,  a modern  study  of  electron  trans- 
port properties  in  HC1  discharges  is  necessary  to  develop  a cross  section  data 
base  having  high  confidence  limits.  For  the  present  analysis,  the  results  of 
Bailey  and  Duncanson  will  be  assumed  to  be  correct;  however,  no  effort  will 
be  made  to  provide  "exact"  agreement  between  the  predicted  and  measured 
transport  properties. 

The  first  reported  effort  directed  towards  developing  a set  of  collision 
cross  sections  for  HC1  which  was  compatible  with  the  available  transport 
data^  was  the  work  of  Nighan  and  co-workers5^  at  UTRC.  This  analysis 
was  performed  using  the  technique  of  Frost  and  Phelps  ^ mentioned  earlier, 
and  included  the  collision  processes  of  momentum  transfer,  rotational 
excitation,  vibrational  excitation,  dissociative  attachment  a nd  electronic 
state  excitation.  At  the  time  of  this  analysis  only  the  cross  section  for 
dissociative  attachment  had  been  measured.  ^ Guidelines  used  in  estimating 
the  remaining  cross  sections  included  the  theoretical  study  of  Itikawa  and 
Takayanagi  ^oon  rotational  excitation  and  elastic  scattering  in  HC1,  and  the 
electron-HCl  excitation  spectrum  reported  by  Compton  et  al.  ^ 

The  cross  section  base  deduced  by  Nighan  et  al  is  shown  in  Fig.  II - 1 
which  is  taken  from  Ref.  54  . The  dominant  feature  of  this  figure  is  the  high 
cross  section,  ~10-^5’  cm^,  for  vibrational  excitation  of  HC1.  It  is  interest- 
ing to  note  that  although  the  onset  energy  for  vibrational  excitation  of  HC1  is 
* 0.36  eV,  the  predicted  onset  for  this  process  from  the  UTRC  analysis  is  at 

~ 0.  46  eV. 

The  availability  of  recent  measurements  ^ on  the  vibrational  ex- 
t citation  cross  sections  for  HC1,  taken  since  the  time  of  the  UTRC  analysis, 

has  provided  motivation  for  a new  determination  of  the  total  cross  section 
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TABLE  II - 1 


MEASURED  TRANSPORT  PROPERTIES  OF  HC  L DISCHARGES 


E/N*10*k  v-cm^ 

* 

SK  ’ SV 

x 10  ^ , cm/sec 

3 

0. 084  - 0.  093 

3.  5 

4.  5 

0.  145 

4.  4 

6 

0.  162  - 0.  172 

4.  5 

9 

0.  294  - 0.  32 

5.  0 

12 

0. 54  - 0.  56 

6.  0 

Lower  values  from  Ref.  51.  The  higher  values  are  reported 
in  Ref.  52  as  the  results  of  Ref.  51.  The  cause  of  this  discrepancy 
is  not  clear. 
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PRELIMINARY  HCI  CROSS  SECTIONS 
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SECTION  - cm 


Fig.  II- 1 


Electron  Collision  Cross  Sections  for  HCI  as  Taken 
F rom  Ref.  54. 
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base  for  HC1  again  by  comparison  with  the  transport  data  of  Bailey  and 
Duncanson.  This  redetermination  has  been  performed  by  a technique 
similar  to  that  of  Frost  and  Phelps,  but  using  a computer  code  for  the 
solution  of  the  Boltzmann  transport  equation  which  was  originally  developed 
by  Carleton  and  Megill^  and  subsequently  modified  by  Lenander^  and  by 
Yos  et  al.1^  The  primary  difference  between  the  two  techniques  is  in  the 
manner  in  which  rotational  excitation  is  evaluated  and  this  is  discussed  in 
Appendix  B. 

Of  the  recent  studies  of  electron  impact  vibrational  excitation  of 

C O _ C Q 

HC1  the  most  detailed  is  that  of  Rohr  and  Linder.  ' Their  preliminary 

C Q > 1 

results  ° for  excitation  of  the  first  vibrational  level  of  HC1,  presented 
only  for  a scattering  angle  of  120°,  exhibited  a sharp  peak  near  threshold 
followed  by  a broad  secondary  peak  centered  around  3 eV.  Although  these 
data  had  not  yet  been  accurately  integrated  over  all  scattering  angles,  the 
cross  section  energy  dependence  did  not  seem  to  vary  significantly  with 
scattering  angle,  and  the  120°  data  was  taken  as  characteristic.  The  angle 
integrated  cross  section  of  the  threshold  peak  was  estimated  to  be 
1.  3 x IQ'15  cm^  ± 50%.  In  actuality,  only  the  relative  cross  section  was 
measured  in  this  experiment.  Their  relative  cross  section  at  an  electron 
energy  of  10  eV  was  normalized  to  an  absolute  measurement  of  the  total 
scattering  cross  section  for  HC1  taken  by  Bruche.^  The  uncertainty  in 
Bruches  measurement  is  not  specified. 

59 

More  recently,  the  differential  cross  section  has  been  integrated 
over  scattering  angle,  and  it  has  been  found  that  the  secondary  peak  takes 
on  more  of  the  appearance  of  a plateau.  Furthermore,  the  magnitude  7 
of  the  cross  section  at  the  threshold  peak  was  raised  to  2.  18  x 10“^  cm  . 
The  uncertainty  in  this  value  is  still  not  defined,  being  primarily  due  to  the 
uncertainty  in  Bruches  measurement.  The  resulting  cross  section  is  shown 
in  Fig.  II - 2 . 

Ziesel  et  al.^  have  also  studied  the  electron  impact  vibration  exci- 
tation of  HC1  near  threshold.  Their  results,  which  span  the  electron  energy 
range  of  0.  36  - 0.  50  eV,  exhibit  a sharp  rise  near  threshold  which  levels 
off  approximately  0.  06  eV  above  onset  in  agreement  with  the  results  of  Rohr 
and  Linder.  In  the  work  of  Ziesel  et  al,  the  vibrational  excita'ion  cross 
section  was  measured  relative  to  that  for  dissociative  attachment  of  HC1. 

As  pointed  out  by  Ziesel  et  al.  , the  three  available  measurements  of  this 
cross  section^’  65,  66  differ  by  a factor  of  five,  with  that  of  Christophorou 
et  al.  ^ ^ being  the  largest.  Thus,  the  peak  cross  section  near  threshold  as 
deduced  from  Ziesel  et  al's.  measurements  can  vary  between  2.  5 x 10"^  - 
1.25  x 10"  ^5  crn2  depending  upon  the  choice  of  the  dissociative  attachment 
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cross  section.  Furthermore,  as  discussed  in  Ref.  60,  an  additional  factor 
of  two  uncertainty  arises  through  the  data  analysis.  Thus,  the  magnitude 
of  the  cross  section  for  vibrational  excitation  of  HC1  cannot  be  readily 
deduced  from  these  data  at  the  present  time.  The  observed  cross  section, 
arbitrarily  normalized  to  a peak  value  of  2.  18  x 10"  cm  , is  shown  in 
Fig.  II  —2  and  can  be  seen  to  be  in  reasonable  agreement  with  the  results 
of  Rohr  and  Linder.  ^9  Although  not  published,  Schultz^  has  independently 
observed  a secondary  peak  at  electron  energies  near  3 eV. 

The  only  other  study  of  the  vibrational  excitation  cross  section  is  the 
work  of  Center  and  Chen,^2  where  the  electron  impact  excitation  rate 
constant  for  HC1  was  determined  in  discharges  of  mixtures  of  HC1-CO-N2- 
Ar.  The  characteristic  energies  in  these  experiments  varied  from  0.  8 - 
1.  7 eV,  and  from  the  measured  rate  constants,  Center  and  Chen^estimated 
that  the  HC1  excitation  cross  section  was  ~ several  x 10'^  cm  . The  data 
taken  in  this  experiment  exhibited  a certain  amount  of  ambiguity  and  was 
re -analyzed  as  part  of  the  present  effort.  This  re -analysis  eliminated  the 
ambiguities,  and  through  detailed  discharge  predictions,  it  was  shown  that 
the  measured  excitation  rate  constants  were  consistent  with  an  HC1  vibra- 
tional cross  section  having  an  energy  dependence  as  measured  by  Rohr  and 
Linder^9  with  a value  near  threshold  of  approximately  10-^  CJ-n  . The 
details  of  this  analysis  are  presented  in  Appendix  C. 

F rom  the  data  discussed  above,  it  may  be  concluded  that  the  energy 
dependence  of  the  electron  impact  vibrational  excitation  cross  section  for 

r 4 

HCl  is  significantly  different  than  that  deduced  in  the  UTRC  analysis, 
(although  this  latter  cross  section  was  consistent  with  the  measured  trans- 
port properties  of  HCl  discharges).  However,  there  is  still  a large  degree 
of  uncertainty  in  the  absolute  magnitude  of  this  excitation  cross  section. 

A new  set  of  trial  cross  sections,  incorporating  these  recent  data  discussed 
above  was  developed  and  used  with  the  previously  described  Boltzmann 
transport  equation  computer  code  to  provide  predicted  transport  properties 
for  HCl.  These  calculations  are  specific  to  a translational  temperature  of 
300°K. 

The  v = 0 -»  1 vibrational  excitation  cross  section  was  taken  to  ha  /e 
the  energy  dependence  as  prescribed  by  Rohr  and  Lindem  1 except  for  he 
near  threshold  region  where  the  results  of  Ziesel  et  al.D®  were  employed. 
This  was  done  because  the  latter  measurements  were  performed  with  finer 


The  situation  is  further  complicated  by  the  fact  that  the  graphical  and 
tabular  data  presented  in  Ref.  60  are  inconsistent  by  a factor  of  two. 
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energy  resolution  than  those  of  Rohr  and  Linder.  There  is  very  little 
information  available  on  the  v = 0 -»  2 vibrational  excitation  cross  section. 
Ziesel  et  al.  report  that  onset  of  this  process  occurs  near  threshold, 
and  that  the  peak  cross  section  is  approximately  one  tenth  that  of  v = 0 -»  1 
excitation.  Rohr  and  Linder  report  that  at  an  energy  of  3 eV  the  ratio  of 
v = 0-*2tov  = 0-»l  cross  sections  is  approximately  one  quarter.  The 
cross  section  used  in  the  present  calculations  was  chosen  to  satisfy  these 
criteria  and  to  have  an  energy  dependence  similar  to  that  observed  for  the 
v = 0 -*  1 cross  section.  The  momentum  transfer  cross  section  used  was 
similar  to  that  developed  in  the  UTRC  analysis^  except  that  the  resonance 
contribution  was  shifted  to  lower  energies  in  concert  with  the  vibrational 
excitation  cross  section.  Indeed,  sharp  structure  in  the  differential  cross 
section  for  elastic  scattering  of  electrons  from  HCl  has  been  observed  to 
occur  near  0.  36  eV,  the  threshold  for  vibrational  excitation.  The  data 
used  for  the  remaining  inelastic  processes  were  the  dissociative  attachment 
cross  section  of  Christophorou  et  al.  ,55  the  electronic  state  excitation 
cross  section  determined  in  the  UTRC  analysis, ^ and  the  ionization  cross 
section  measured  by  Compton  and  Van  Voorhis.^  Rotational  excitation 
was  handled  by  the  single  level  approximation  as  described  in  Appendix  B. 


The  above  set  of  cross  sections  was  adjusted  until  "reasonable" 
agreement  was  obtained  between  measured  and  predicted  transport  proper- 
ties. For  the  most  part  these  adjustments  were  limited  to  variations  in 
the  magnitude  of  the  vibrational  excitation  cross  sections  and  to  the  shape 
and  magnitude  of  the  momentum  transfer  cross  section  in  the  resonance 
region.  The  final  "best  set"  of  vibrational  and  momentum  transfer  cross 
sections  is  shown  in  Fig.  H-3.  Note  that  the  peak  value  of  the  vibrational 
excitation  cross  section  is  1.  3 x 10-13  cm  which  is  in  the  midrange  of  the 
available  estimates . 58-62  Predicted  transport  properties  for  HCl, 

developed  from  the  cross  sections  shown  in  Fig.  II-3,  are  compared  with 
the  measurements  of  Bailey  and  Duncanson^  in  Tables  H-2  and  H-3. 

Shown  for  comparison  are  the  predictions  from  Ref.  54.  As  can  be  seen 
the  present  calculated  results  are  reasonable,  although  the  predicted  values 
of  and  Vj)  are  somewhat  higher  than  the  data  at  the  highest  values  of 
E/N  considered.  Although  much  better  comparison  could  have  been  achieved 
by  further  adjustment  of  the  cross  sections,  the  accuracy  of  the  transport 
data  does  not  warrant  the  effort.  Interestingly  enough,  the  difference 
between  the  present  predictions  and  those  of  Ref.  54  are  not  large  although 
the  cross  sections  used  are  quite  different. 


Calculations  have  also  been  performed  for  mixtures  of  10%  HCl/90% 
Ar  and  3.33%  HCl/96.  67%  Ar.  Standard  values  for  the  Ar  cross  sections 
were  used,  as  listed  in  Ref.  4.  The  predicted  characteristic  energies  and 
drift  velocities  for  these  mixtures  are  shown  in  Fig.  II-4  vs  E/N^^. 
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Fig.  II- 


Electron  Collision  Cross  Sections  for  HC1  as  Selected 
for  the  Present  Study. 
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TABLE  H-2 

COMPARISON  OF  PREDICTED  AND  MEASURED  CHARACTERISTIC 
ENERGIES  IN  HCL  DISCHARGES 


eK’ 

eV 

Present 

Analysis 

E/N  x 10+*  ,V  -cm^ 

Data,  Ref.  51 

Prediction,  Ref.  54 

3 

0.  084-0.  093 

0.09 

0.  08 

4.  5 

0.  145 

, * 

(0. 16) 

0.  12 

6 

0.  162-0.  175 

(0.21) 

0.  165 

7 

— 

0.24 

0.20 

9 

0.294-0.  32 

(0.36) 

0.  33 

10 

— 

0.43 

0.43 

12 

0.54  - 0.  56 

(0.74) 

0.  73 

15 

— 

1.32 

1.  31 

( ) -»  Interpolated  from  Predictions  Presented  in  Ref.  54 
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TABLE  H-3 


:V 


COMPARISON  OF  PREDICTED  AND  MEASURED  DRIFT 
VELOCITIES  IN  HCL  DISCHARGES 


VD*  106 

, cm/sec 

E/Nx  10+l6,V-cm2 

Data,  Ref.  51 

Prediction,  Ref.  54 

3 

2.  5 

2.  7 

4.  5 

4.4 

(4.4)* 

6 

4.  5 

(4.  95) 

7 

— 

5.4 

9 

5.0 

(6.0) 

10 

— 

6.  1 

12 

. 6.  0 

(6.25) 

15 

— 

6.  4 

Present  Analysis 

2.5 

3.  7 

4.  5 

5.  1 
5.7 

5.  8 

6.  1 

6.4 


( ) -*  Interpolated  from  Predictions  Presented  in  Ref.  54 


/sec 


The  results  for  pure  HC1  are  also  shown  for  purposes  of  comparison. 

As  can  be  seen,  the  presence  of  Ar  affects  the  transport  properties  at  the 
higher  values  of  where  both  the  drift  velocity  and  characteristic 

energy  decrease  with  increasing  Ar  concentration.  This  phenomenon  is 
primarily  due  to  momentum  transfer  collisions  with  Ar  and  is  less 
pronounced  at  the  lower  values  of  E/N  (and  thus  lower  values  of  ) 
because  of  the  deep  Ramsauer  minimum  occuring  in  the  Ar  momentum 
transfer  cross  section  at  low  electron  energies. 

The  most  interesting  feature  of  the  discharge  predictions  from  the 
viewpoint  of  laser  performance,  is  the  manner  in  which  the  discharge 
energy  is  distributed  among  the  various  inelastic  processes  occuring 
in  the  gas.  The  percentage  energy  allocation  into  rotational,  vibrational 
and  electronic  state  excitations  as  a function  of  is  shown  in 

Fig.  H-5  for  the  cases  of  pure  HC1  and  3.33 % HC1/96.  67%  Ar.  The 
predictions  for  the  10%  HCl/90%  Ar  case,  not  shown,  are  similar  to  those 
for  the  3.33%  mixture.  As  can  be  seen  for  E/Npid  > 8 x 10"^  V-cm^, 
approximately  80  - 90%  of  the  discharge  energy  is  channeled  into  vibra- 
tional excitation.  Furthermore,  the  addition  of  Ar  is  beneficial  in  that  it 
extends  the  range  of  Townsend  parameter  over  which  a high  level  of  vibra- 
tional excitation  may  be  achieved,  while  simultaneously  decreasing  un- 
desirable losses  to  electronic  state  excitation. 

The  actual  rate  constants  for  electron  impact  excitation  of  the  first 
and  second  vibrational  levels  of  HC1,  dissociative  attachment  and  electronic 
state  excitation  are  shown  vs  E/N  in  Figs.  II-6  to  8 for  the  cases  of  pure 
HC1,  10%  HCl/90%  Ar  and  3.33%  HCl/96.67%  Ar,  respectively.  Note  that 
the  vibrational  excitation  rate  constant  is  of  order  10-®  cc/sec,  and  that 
at  the  higher  values  of  E/N,  the  rate  constant  for  direct  excitation  to  v = 2 
is  as  much  as  15%  of  that  for  excitation  to  v = 1.  Furthermore,  it  can  be 
seen  that  the  rate  constant  for  dissociative  attachment  peaks  in  the  same 
region  as  that  for  vibrational  excitation  and  has  a value  of  ~10"^®  cc/sec. 

It  should  be  remembered  that  the  magnitude  of  the  cross  section  for  this 
process  is  still  controversial,  and  that  the  largest  measured  value  ^ has 
been  used  in  the  present  analysis.  Lastly,  it  can  be  seen  from  comparison 
of  Figs.  U-6  and  II-7  that  the  presence  of  Ar  severely  degrades  electronic 
state  excitation  without  significantly  altering  the  rate  constants  for  vibra- 
tional excitation  and  dissociative  attachment.  This  is  consistent  with  the 
results  presented  in  Fig.  II - 5 . 
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Fig.  II-7  Predicted  Electron  Excitation  Rate  Constants  for  HCl  in 
a 10%  HCl/90%  Ar  Mixture. 
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Fig.  II- 8 Predicted  Electron  Excitation  Rate  Constants  for  HC1  in 
a 3.  33%  H Cl/96.  67%  Ar  Mixture. 
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c. 


LASER  PERFORMANCE  STUDIES 


1 . Gain  Predictions 

The  results  of  the  previous  section  have  been  coupled  with  a vibra- 
tional kinetics  code  in  order  to  provide  gain  predictions  for  generic 
experimental  conditions  studied  in  the  NRL  program.  The  vibrational 
kinetics  package  employed,  discussed  in  Section  I.C.  1 of  this  report, 
included  all  vibration-vibration  and  vibration-translation  reactions  of 
importance  in  HCl/Ar  mixtures.  The  calculations  were  Performed  under 
the  assumption  of  a uniform  medium  having  an  active  lergth  of  100  cm. 
Rotational  equilibration  with  translation  was  assumed  to  be  instantaneous, 
so  that  the  energy  supplied  to  the  gas  by  electron  impact  rotational  excitation 
was  equi -partitioned  between  the  translational  and  rotational  degrees  of 
freedom. 

A specific  case  considered  corresponded  to  a mixture  of  10% 

HCl/90%  Ar  at  a pressure  of  one  atmosphere  and  initial  temperature  of 
300°K.  An  E/N  of  1.2  x 10  V-cm^  was  selected,  a value  corresponding 
to  a predicted  high  efficiency  for  vibrational  excitation  (see  Fig.  II _ 5 ) . 

The  sustainer  current  was  taken  to  be  6 amps/cm^,  a value  typical  of 
those  realized  experimentally.  With  a predicted  drift  velocity  of 
5.  5 x 106  cm/sec  this  sustainer  current  would  correspond  to  an  electron 
density  of  ~7  x 10  /cm^.  The  two  vibrational  excitation  processes  in- 

cluded in  the  calculation  were 


e + HC1  (v)  -»  e + HC1  (v  + 1) 

a~d  (II.  1) 

e + HC1  (v)  -»  e + HC1  (v  + 2) 

- 8 - 9 3 

having  rate  constants  of  1.3  x 10“  and  1.3  x 10  cm  /sec,  respectively. 
The  e-beam  pulse  was  taken  to  have  a duration  of  3 psec. 

The  predicted  small  signal  gain  histories  for  this  case  are  shown 
in  Fig.  II-9  for  the  lowest  vibrational  levels  of  HC1.  As  can  be  seen,  the 
predictions  exhibit  small  signal  gains  greater  than  10-^  cm-1.  The  de- 
creasing gain  shortly  after  pulse  termination  is  a result  of  vibrational- 
translational  deactivation  with  a commensurate  increase  in  translational 
temperature.  The  total  predicted  increase  in  translational  temperature 
over  20  psec  was  ~ 80°K. 
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Calculation  of  Small  Signal  Gain  of  Several  HC1  Vibrational 
Transitions  as  a Function  of  Time  After  the  Initiation  of  an 
E-Beam  Sustainer  Discharge  in  a 10%  HCl/90%  Ar  Mixture. 
The  Assumed  Experimental  Conditions  are  Listed  on  the 
Figure.  No  H or  Cl-atoms  are  Assumed  to  be  Present. 


Unfortunately,  such  high  gains  were  not  realized  experimentally. 
Indeed,  as  discussed  in  detail  in  Ref.  69,  although  a large  range  of  mixture 
ratios,  pressures  and  E/N's  were  studied,  lasing  was  never  observed  to 
occur,  suggesting  that  actual  small  signal  gains  were  < 10"3  cm"*.  This 
observation  led  to  a review  of  the  assumptions  used  in  the  modeling. 

One  feature  that  had  not  been  considered  was  the  possible  importance  of 
atomic  species,  created  in  the  discharge,  which  could  vibrationally  de- 
activate HC1. 

Atomic  species  are  created  directly  in  the  discharge  via  dissociative 
attachment,  i.  e.  , 

« 

e + HCl  -»  H + Cl"  . (H -2) 

Under  optimum  operating  conditions  for  vibrational  excitation,  the  discharge 
is  attachment  dominated, and  thus  process  (II-2)  leads  to  the  dissociation  of 
one  HCl  molecule  for  every  secondary  electron  created.  An  additional 
source  of  atomic  species  in  the  discharge  is  the  reaction 

Ar+  4 HCl  ->  Ar  H+  + Cl  , (II -3) 

which  has  recently  been  reported70  to  have  a rate  constant  > 1.  3 x 10-1° 
cm3/sec.  This  reaction  acts  in  competition  to  the  Ar+  clustering  reaction 


Ar+  + 2 Ar  -*  Ar^+  + Ar  (II -4) 

T 1 316  1 

which  has  a rate  constant  of  ~ 10  cm  -sec  . The  reaction  between 
Ar£  and  HCl  has  not  been  studied,  however,  the  branch  to  form  ArH+ 
would  be  close  to  thermoneutral  and  is  perhaps  not  favorable.  For  a 10% 
HCl/90%  Ar  mixture  at  1 atmosphere  reaction  (II- 3)  is  favored  over 
(II-4),  and,  thus,  approximately  one  chlorine  atom  is  created  for  every 
positive  ion.  In  the  instance  when  ArH+  is  the  dominant  ion,  it  most 
probably  reacts  with  Cl"  to  reform  HCl,  i.  e. , 

ArH1  + Cl"  ->  HCl  + Ar  , (II -5) 

and  no  additional  free  atom  would  be  formed.  Under  conditions  where  Ar ^ 
is  the  dominant  ion,  it  could  recombine  with  Cl"  to  form  free  Cl.  Neither 
of  these  last  two  reactions  have  been  studied. 
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Thus,  it  would  appear  that  in  the  attachment  dominated  discharge, 
at  least  two  atoms  will  be  formed  for  every  ion  pair.  For  the  beam 
operating  conditions  of  V = 300  kV  and  I m 6 amps /cm*3,  the  ion  pair 
production  rate  is  predicted  to  be  ~ 3 x lO^Vcm^-sec.and  thus  for  a 
3 psec  pulse  > 1.8x  10^,  atoms  will  be  formed.  This  is  to  be  compared 
with  an  HC1  concentration  of  2.  5 x 10l8/crn3  for  a 10%  HC1  mixture. 

From  these  two  concentrations  it  can  be  seen  that,  if  the  rate  constant 
for  vibrational  deactivation  of  HC1  by  atoms  is  >150  times  that  for  self 
deactivation,  the  atoms  will  provide  the  dominant  vibrational  loss 
mechanism. 

72 

The  rate  constants  for  vibrational  deactivation  of  HC1  (v  = 1)  by  H 
and  Cl73  have  been  reported  to  be  6.  5 and  8.  5 x 10“^  cm3 /sec,  respec- 
tively, at  room  temperature.  This  is  to  be  compared  to  a rate  constant  of 
2.4  x 10“*^  cm^/sec  for  self  deactivation  of  HC1,  again  at  room  temperature. 
Thus,  a concentration  of  1.8  x 10^°  H and  Cl  atoms  is  sufficient  to  more 
than  triple  the  overall  HC1  vibrational  deactivation  rate  in  a 1 atm  mixture 
of  10%  HCl/90%  Ar.  Gain  predictions  have  been  performed  for  the  identical 
conditions  of  Fig.  II -9  but  with  three  times  the  HC1  deactivation  rate 
constant.  The  results  of  this  calculation  are  shown  in  Fig.  11-10,  and  it 
can  be  seen  that  the  increased  deactivation  severely  degrades  the  small 
signal  gain  to  levels  which  would  be  only  marginally  observable  in  the 
experiment. 

Thus  it  would  seem  that  the  presence  of  discharge  created  H and  Cl 
atoms  in  the  gas  can  inhibit  laser  action.  There  are  several  other  possible 
sources  of  atom  production  which  have  not  been  considered  in  the  above 
analysis.  For  example,  in  the  case  of  the  primary  beam,  HC1  can  be 
directly  ionized,and  the  HC1+  can  later  recombine  dissociatively. 
Furthermore,  direct  dissociation  or  dissociative  ionization  can  occur. 

For  the  secondary  electrons,  electronic  state  excitation  leading  to  dissocia- 
tion can  be  important  at  higher  values  of  E/N  as  demonstrated  in  the 
discharge  calculations  presented  earlier.  Lastly,  any  Ar  metastable  states 
produced  in  the  discharge  could  be  collisionly  quenched  by  HC1  in  a dissocia- 
tive process.  Such  processes  can  only  be  conjectured  at  present  because 
of  the  limited  experimental  evidence  available. 

2 . Sustainer  Current  Measurements 

If  a high  concentration  of  atomic  species  limits  laser  performance, 
one  could  envision  either  operating  under  conditions  which  minimize  their 
importance  or,  alternately,  introducing  atomic  scavengers  into  the  gas 
mixture.  Unfortunately,  there  is  another  aspect  of  the  experimental 
measurements  which  is  not  understood  at  present;  this  has  to  do  with  the 
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Fig.  II- 10  Calculation  of  Small  Signal  Gain  of  Several  HCl  Vibrational 
Transition  as  a Function  of  Time  After  the  Initiation  of  an 
E-Beam  Sustainer  Discharge  in  a 10%  HCl/90%  Ar  Mixture. 
The  Experimental  Conditions  are  Assumed  Identical  to  those 
of  Fig.  II— 9 Except  that  to  Simulate  the  Presence  of  Atoms 
Produced  in  the  Discharge,  the  HCl  V-*T  Rate  has  been 
Multiplied  by  3. 
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the  observed  sustainer  currents.  These  currents  are  not  only  lower  than 
theoretically  predicted  but  also  behave  anomolously  with  varying  E/N, 
pressure  and  mixture  ratio.  As  a specific  example,  in  an  attachment 
dominated  discharge  the  local  electron  density  is  specified  by  the  relation- 
ship 


N = 
e 


a 


k Na 
a A 


(II -6) 


where  a is  the  ion  pair  production  rate,  ka  is  the  attachment  rate  constant 
and  is  the  number  density  of  the  attaching  gas.  For  10%  HCl/90%  Ar 
at  1 atm  with  an  a of  3 x 1 0^  1 /cm-* -sec,  characteristic  of  the  NRL  e-beam 
operating  conditions,  Eq.  (II -6)  reduces  to 


N 

e 


1.2  x 103/k  . 

a 


(II -7) 


This  electron  density  corresponds  to  a sustainer  current  defined  by 


is  = e Ne  VD 


= 1.9  x 10"16  V^/k  , 
u a 


(II -8) 


and,  thus,  the  sustainer  current  may  be  predicted  from  the  transport 
properties  given  in  Figs.  U-4  and  II -6 . 

As  discussed  in  detail  in  Ref.  69,  the  measured  sustainer  current 
was  always  lower  than  that  predicted  from  Eq.  (II-8).  A typical  comparison 
is  shown  in  Fig.  11-11  which  was  excerpted  from  Ref.  69.  Measurements 
were  extended  to  other  pressures  and  mixture  ratios  with  equally  diverse 
results.  For  example,  if  the  value  of  E/N  and  the  mixture  ratio  is  held 
fixed  and  the  pressure  increased,  the  electron  density  and  jg  should  remain 
constant.  This  is  because  both  a and  Na  scale  linearly  with  pressure. 

The  experimental  observation  was  that  j8  decreased  with  increasing 
pressure. 

At  first  it  was  felt  that  this  discrepancy  between  theory  and  experi- 
ment could  have  been  a result  of  additional  attachment  processes  occuring 
in  HC1.  For  example,  apparent  attachment  effects  have  been  observed 
near  thermal  energy  in  HC1  although  they  are  not  properly  understood 
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Fig.  II- 11  Comparison  of  Predicted  and  Observed  Sustainer  Current 
for  an  E-Beam  Sustainer  Discharge  in  a 10%  HCl/907o  Ar 
Mixture  at  a Total  Pressure  of  1 Atmosphere. 


(e.g.,  see  Ref,  55  and  the  discussion  in  Ref.  74).  To  examine  this 
possibility, sustainer  current  measurements  were  performed  in  pure  Ar 
and  N2.  Again  unexpected  behavior  was  observed. ^9  The  pressure 
dependence  of  the  measured  sustainer  current  in  N£  is  shown  in  Fig.  11-12 
for  several  values  of  E/N.  Under  the  experimental  conditions  the  N2  dis- 
charge is  expected  to  be  controlled  by  dissociative  recombination^^'  7^ 
with  the  dominant  ion  being  . In  this  latter  case  it  can  be  shown^’  75 
that  at  constant  E/N  the  sustainer  current  should  scale  as  the  square  root 
of  pressure.  Indeed  such  a pressure  dependence  has  been  observed  in  N2 
between  0.2-1  atm  in  another  e-beam  sustainer  discharge  experi-nent, 
although  under  different  beam  conditions  than  used  in  the  present  study. 

As  can  be  seen  from  Fig.  11-12,  a much  more  complicated  pressure 
dependence  was  observed  in  the  NRL  experiment. 

Thus,  the  observed  behavior  of  sustainer  current  is  not  understood 
at  present.  Particularly  in  the  case  of  N2,  it  is  difficult  to  conceive  of  a 
reasonable  kinetic  electron  loss  mechanism  which  could  produce  the 
observed  pressure  dependence.  Most  probably  these  complications  are  in 
the  ion  source  term,  caused  perhaps  by  beam  spread  and  non-uniform 
energy  deposition.  Additional  diagnostics,  not  presently  available  on  the 
experimental  facility,  would  be  required  to  examine  this  possibility.  The 
possible  effects  of  trace  impurities  in  the  test  gases  and  medium  in- 
homogeneity must  also  be  cons ide red. ^ 

3.  Energetics 


Before  closing  the  discussion  on  laser  performance,  there  is  an 
additional,  somewhat  subtle  point  which  deserves  mention.  In  general,  in 
the  case  of  e-beam  sustained  infrared  lasers,  the  amount  of  energy 
deposited  into  the  gas  by  the  electron  beam  is  negligible  in  comparison  to 
the  energy  deposited  by  the  sustainer  field.  This  is  not  the  case  in  the 
present  experiment.  A 300  kV  electron  penetrating  an  atmosphere  of  Ar 
loses  approximately  2.8  kV/cm.7^  For  an  electron  beam  current  of 
6 amps/cm^  this  corresponds  to  an  energy  deposition  rate  of  17  k watts /cm3. 
In  contrast, from  Fig.  II -6  it  can  be  seen  that  in  a 10%  HCl/90%  Ar  mixture 
at  one  atmosphere  a sustainer  current  of  ~ 7 amps /cm^  is  attained  with  a 
sustainer  field  of  2.  3 kV/cm,  corresponding  to  a sustainer  energy  input 
rate  of  16  k watts  /cm  . For  this  case,  beam  and  sustainer  energy  inputs 
are  approximately  equal.  Although  the  sustainer  energy  is  primarily  con- 
verted to  vibrational  excitation,  the  beam  energy  degrades  predominantly 
to  translational/rotational  heat  content.  Indeed  for  a 3 ns  e-beam  pulse  an 
energy  input  rate  of  16  kW/cm^  is  sufficient  to  raise  the  gas  temperature 
80°K . The  gas  heating  caused  by  electron  beam  energy  deposition  was  not 
included  in  the  gain  predictions  presented  in  Section  II.  C.  1.  The  effect  of 
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Fig.  11-12  Measured  Sustainer  Currents  as  a Function  of  Pressure  for 
Various  Values  of  E/N  in  Pure  N^.  In  the  Lower  Right-Hand- 
Side  of  the  Figure  are  Lines  Indicating  the  Expected  Slope  for 


the  increased  gas  temperature  would  be  to  decrease  the  small  signal  gain 
predictions  below  the  levels  shown  in  Fig.  II-IO. 

More  importantly,  inasmuch  as  only  a small  fraction  of  the  beam 
power  is  deposited  in  the  gas,  the  low  sustainer  powers  would  suggest  that 
at  best  the  HC1  EDL  would  be  very  inefficient.  For  example  the  total 
beam  energy  flux  at  300  kV  and  6 amps /cm2  is  1. 8 M watts /cm2.  In 
comparison  a sustainer  field  of  2.  3 kV/cm  applied  over  10  cm  with  a 
sustainer  current  of  7 amps /cm2  corresponds  to  a sustainer  energy  flux 
of  0.  16  MW /cm2.  Even  if  all  the  sustainer  energy  could  be  converted  to 
laser  photons,  the  laser  would  be  less  than  10%  efficient.  The  primary 
way  to  improve  the  maximum  efficiency  is  to  increase  the  sustainer  current. 
If  the  discharge  were  truly  uniform  and  attachment  dominated,  this  could 
be  achieved  by  decreasing  the  attachment  rate  either  by  decreasing  the 
partial  pressure  of  HC1  or  by  operating  at  E/N's  where  the  attachment  rate 
constants  are  predicted  to  be  low.  Alternatively,  one  could  increase  the 
sustainer  current  by  increasing  the  total  gas  pressure  while  keeping  the 
HC1  partial  pressure  constant.  Unfortunately,  as  discussed  in  Section  II.  C.  2 
and  in  Ref.  69,  the  sustainer  current  did  not  behave  in  the  expected  manner 
when  such  variations  were  examined  experimentally. 


-106- 


D.  INFRARED  RADIATION  STUDIES 

At  the  beginning  of  the  HC1  EDL  program  attempts  were  made  to 
deduce  information  about  the  discharge  and  vibrational  kinetics  by  monitor- 
ing the  temporal  fluorescence  histories  of  the  vibrational/rotational  bands 
of  HC1.  This  study,  which  was  not  successful,  has  been  described  in  Ref. 77 
and  will  only  be  discussed  briefly  here. 

The  first  measurements  were  made  in  the  3-4 p wavelength  region 
utilizing  an  indium  antimonide  detector  so  positioned  as  to  observe  radiation 
emanating  from  the  full  100  cm  active  length  of  the  device.  Initially,  funda- 
mental band  fluorescence  measurements  were  performed  in  discharges  of 
10%  HCl/90%  Ar.  However,  the  gas  was  found  to  be  severely  optically  thick 
under  such  conditions,  and  the  measurements  were  not  useful. 

The  problem  of  optical  thickness  in  HCl/Ar  mixtures  was  examined 
in  some  detail.  It  can  be  shown^  for  an  isolated  Lorentz  broadened 
spectral  line,  that  the  fraction  of  total  emitted  radiation  which  will  escape 
through  a radiating /absorbing  gas  of  length  L is  given  by 


I/I  4--  ii  = F <x)/x  » 

optically  thin 


(II  -9 ) 


where 


x = S [HC1]  L/tt  Av  t = k , L/2 

1j  C • 1 • 


(II  - 10) 


- 1 -2 

S is  the  integrated  line  strength  in  ama  -cm  , [HCl]  is  the  concentration 
of  HCl  in  amagats,  Avl  is  the  full  Lorentz  broadened  line  width  at  half 
height  in  cm~I  and  kc  is  the  absorption  coefficient  at  line  center.  The 
function  F(x)  is  defined  by  the  relationship 


F (x ) = x e ~X  [ JQ  (ix ) - i J 1 (x ) ] , (II- 

where  the  J's  are  Bessel  functions.  This  expression  simplifies  somewhat 
for  large  values  of  x,  where  it  can  be  shown  78  that 


F(x)  » (2x/tt) 


1 


I 

t 

p 
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and  thus  that 

1 ^optically  thin  **  (ttx)  ( TTk  L ) (11-13) 

' 7 ' c.  1.  7 
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Typical  values  of  the  R and  P branch  v = l-»0  line  strengths  of 
HC1  are  listed  in  Table  II-4  along  withArgon  broadened  line  widths®0  and 
the  evaluated  centerline  absorption  coefficient.  From  Eq.  (II - 1 0 ) and 
Table  II -4,  it  can  be  seen  that  for  a mixture  of  10%  HCl/90%  Ar  at  1 atm 
and  L = 100  cm  the  quantity  x is  quite  large,  and,  thus,  the  approximate 
expression  (11-13)  may  be  used.  The  fractionally  transmitted  radiation  for 
these  conditions  is  ~0.04  for  significantly  populated  rotational  lines. 

The  full  expression  (II-9),  has  been  evaluated  as  a function  of  the 
quantity  [HCl]  L for  a centerline  absorptio n coefficient  of  89  cm"',  a 
value  corresponding  to  the  strongest  P branch  transitions  of  HCl.  The 
results  are  shown  in  Fig.  11-13,  and  it  can  be  seen  that  conditions  corres- 
ponding to  [HCl]  L<  10-1  are  required  to  minimize  the  effects  of  self 
absorption.  In  actuality  the  predictions  shown  assume  that  the  chlorine  is 
all  isotope  35.  Actually  the  isotope  35  to  37  ratio  is  3 to  1,  and,  thus,  in 
an  optically  thick  environment  radiation  from  HCl®7  will  be  enhanced  over 
that  from  HCl.  3 5 

Since  the  experimental  facility  could  not  be  readily  modified  to 
provide  fluorescence  measurements  over  a shorter  optical  path,  discharge 
experiments  were  performed  in  .01%  - 1%  HCl/Ar  mixtures.  The 
fluorescence  observed  in  these  experiments  exhibited  a decay  rate  (upon 
termination  of  the  e-beam)  orders  of  magnitude  more  rapid  than  would  be 
expected  for  HCl  vibrational  deactivation.  Additional  measurements, 
performed  with  the  emission  screened  with  an  HCl  gas  filter,  confirmed 
that  the  observed  radiation  was  a result  of  Ar  rather  than  HCl  fluorescence. 

Additional  isea surements,  in  both  pure  Ar  and  Ar/HCl  mixtures, 
were  performed  over  the  3 - 4 y region  using  a series  of  narrow  bandpass 
filters  in  order  to  determine  the  spectral  variation  of  the  radiation.  It 
was  found  that  the  radiation  was  strongly  peaked  at  3.6  and  4.  0 pm,  wave- 
lengths corresponding  to  the  position  of  previously  observed  Ar  lines. 

As  a last  attempt  to  monitor  HCl  fluorescence,  a series  of  measurements 
were  pe rformed  at  wavelengths  of  1.5-  2.  5 pm  in  the  region  of  the  HCl 
first  overtone.  Again,  the  observed  fluorescence  did  not  exhibit  behavior 
characteristic  of  HCl.  A large  number  of  Ar  lines  have  been  reported82 
to  occur  between  1.5  - 2.  5 pm,  and  these  are  most  probably  the  source 
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of  the  observed  radiation.  At  this  point  this  aspect  of  the  experimental 
program  was  abandoned. 

It  is  of  some  interest  to  speculate  on  the  origin  of  the  Argon  fluores- 
cence. The  transitions  observed  occur  between  highly  excited  Ar  states 
which  can  be  excited  directly  by  the  electron  beam.  However,  the 
fluorescence  decay  was  found  to  occur  over  times  as  long  as  5 - 10  usee 
after  termination  of  the  e-beam.  This  is  the  same  order  time  scale  as  for 
electron-ion  recombination  for  the  experimental  conditions.  The  dominant 
positive  ion  in  the  pure  Ar  and  very  dilute  HCl/Ar  mixtures  is  ArJ  and 
thus  the  observed  fluorescence  may  have  been  a result  of  the  dissociative 
recombination  reaction 


+ * 

Ar^  + e -»  Ar  + Ar 


(11-14) 


82 

which  is  known  to  produce  electronically  excited  Ar.  A more  detailed 
investigation  would  be  required  in  order  to  confirm  this  latter  possibility. 


FUNDAMENTAL  BAND 


Fig.  11-13  Relative  Intensity  of  HC1  Fundamental  Vibration/Rotation 
Band  for  Various  Values  of  Optical  Depth.  The  Lines  are 
Assumed  to  be  Pressure  Broadened. 


/[HC1]  = Centerline  absorption  coefficient  for  Ar  broadening; 
normalized  by  the  HC1  mole  fraction 


E.  SUMMARY  AND  CONCLUSIONS 


In  summary,  the  HC1  EDL  modeling  described  above  has  produced 
several  salient  but  contradictory  facts  regarding  possible  system  perfor- 
mance. In  particular  the  HCl/Ar  discharge  predictions  suggest  that  some 
80  - 90%  of  the  discharge  energy  may  be  channeled  into  vibrational  excitation 
of  the  HC1.  Unfortunately,  however,  the  subsequent  kinetic  analysis  has 
demonstrated  that  atomic  species  produced  in  the  discharge,  primarily 
by  dissociative  electron  attachment,  can  severely  degrade  laser  perfor- 
mance through  vibration-translation  de-activation  of  the  HC1.  Since 
dissociative  attachment  and  vibrational  excitation  of  HC1  occur  in  approx- 
imately the  same  electron  energy  range,  the  production  of  atomic  species 
in  the  discharge  cannot  be  completely  eliminated;  however,  it  can  be 
minimized.  A modern  study  of  the  transport  properties  and  attachment 
phenomena  occuring  in  HC1  electron  swarms  would  be  very  useful  in 
determining  optimum  discharge  conditions  for  laser  performance.  The 
use  of  chemical  additives  as  scavengers  for  atomic  species  should  also 
be  considered. 

An  additional  feature  peculiar  to  the  atmospheric  HC1  EDL  is  the 
projected  low  efficiency  of  the  device.  As  discussed  in  the  main  text,  the 
unexpected  variation  of  sustainer  current  as  a function  of  pressure  and 
mixture  ratio  precluded  any  specification  of  the  operational  parameters 
required  for  optimum  conversion  of  electrical  energy  to  vibration.  It  is 
possible  that  the  sustainer  current  observations  were  device  specific. 

Similar  studies  performed  on  a small,  low  current  device,  which  would 
allow  both  absorption  free  fluorescence  measurements  and  local  current 
density  determinations,  would  be  invaluable  in  specifying  the  true  potential 
of  an  HC1  EDL. 
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APPENDIX  A 


Electron  Impact 


Rotational  Excitation  of 


H. 


Homonuclear  molecules  having  non-zero  nuclear  spin  may  have 
both  even  and  odd  rotational  levels  corresponding  to  either  symmetric 
or  antisymmetric  alignments  of  nuclear  spin  respectively,  depending  on 
the  ground  electronic  state  configuration  of  the  molecule.  Since  transitions 
between  symmetric  and  anti- symmetric  alignments  are  strongly  forbidden, 
the  total  concentration  of  odd  and  of  even  states  remains  constant  indepen- 
dent of  temperature  and  is  dependent  only  upon  the  alignment  degeneracies, 
i.  e.  statistical  weights.  The  state,  i.  e.  , odd  or  even,  which  has  the 
higher  statistical  weight  is  termed  the  ortho  modification,  while  the  state 
with  lower  weight  is  termed  the  para  modification. 

In  the  case  of  H , which  has  nuclear  spin  1/2,  the  statistical 
weight  of  odd  states  is  3 and  that  of  even  states  is  1.  Thus,  at  all  tempera- 
tures 75%  of  the  H2  molecules  will  be  in  odd  rotational  levels  and  2 5%  in 
even.  The  situation  is  reversed  in  D^,  which  has  nuclear  spin  1.  For  this 
latter  situation,  even  states  have  statistical  weight  6 and  odd  states  3. 

Thus  66.  7%  of  the  D2  molecules  will  be  in  even  rotational  levels  and  33.  3% 
in  odd.  The  resulting  predicted  rotational  populations  of  the  lower  levels 
of  H2  and  D2  are  shown  vs  temperature  in  Figs.  A-l  and  A -2.  The 
rotational  constants  used  for  these  calculations  were  taken  from  Herzberg.  * 
As  can  be  seen,  at  low  temperatures,  ~ 100°K,  only  the  J = 0 and  1 states 
have  significant  populations,  although  some  18%  of  the  D2  molecules  are  in 
the  J = 2 state.  However,higher  rotational  levels  begin  to  play  a more 
important  role  with  increasing  temperature.  Indeed.by  300°K  4 or  5 
rotational  levels  have  populations  > 10%.  As  can  be  seen  from  Fig.  A -2, 
the  temperature  variation  is  particularly  pronounced  in  the  case  of  D2;  at 
a temperature  of  300°K  J = 2 is  the  most  populated  state  and  the  populations 
of  J = 3 and  4 are  within  a factor  of  two  of  those  for  J = 0 and  1. 

Based  on  these  significant  population  variations  it  might  be  expected 
that  the  amount  of  discharge  energy  absorbed  by  rotation  and  also  the 
transport  properties  of  H2/D2  discharges  would  vary  with  temperature, 
since  each  rotational  state  has  a different  onset  energy  and  excitation  cross 
section.  Crompton^  et  al,  performed  a careful  study  of  the  transport 
properties  of  normal  and  parahydrogen  at  77°K  and  normal  hydrogen  at 
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300°K.  At  77°K  parahydrogen  is  almost  completely  in  the  J - 0 state, 
whereas  normal  hydrogen  is  populated  75%  in  J = 1.  At  300°K,  normal 
hydrogen  has  rotation  populations  of  13,  2,  66.  2,  11.  6 and  8.  7%  in  J = 0 
to  3 respectively.  Thus,  these  three  cases  bracket  a wide  range  of 
rotational  population  distributions.  No  significant  differences  in  the 
transport  properties  of  these  three  cases  were  noted  at  E/N's  corresponding 
to  characteristic  energies  well  above  the  onset  energy  for  rotational 
excitation.  In  particular,  the  observed  differences  were  <>  1%  for  E/N  > 

5 x 10" V-cm  . Crompton  et  al,  ^ point  out  that  this  observation  is  a 
result  of  the  fact  that,  well  above  onset  threshold,  the  rotational  excitation 
cross  sections  scale  inversely  with  the  threshold  energy.  Thus  the 
energy  loss  per  collision  to  rotational  excitation  will  be  independent  of 
initial  rotational  state  for  sufficiently  high  electron  energies. 

This  point  has  been  examined  in  more  detail  in  the  present  analysis. 
Chang  and  Temkin^  have  shown  that  well  above  threshold  the  rotational 
excitation  cross  section  o (J  -*  J1)  scales  ac  the  square  of  the  Clebsch- 
Gordon  coefficient,  C (J2  J':000),  where  J and  J'  are  lower  and  upper 
rotational  states*  respectively, of  the  transition.  For  transitions  J -»  J + 2, 
the  excitation  cross  section  is  proportional  to^ 


^n-cWooo 


(A-  1) 


Furthermore  the  energy  loss  associated  with  each  transition  is 


J + 2 = (4J  + 6,B  ’ 


(A-2) 


where  B is  the  rotational  constant  of  the  molecule.  Combining  Eqs.  (A-l) 
and  ( A-2  ) it  can  be  readily  shown  that  for  both  and  D2  at  electron 
energies  well  above  threshold 


a (J  •*  J + 2)  AEJ-»  J + 2 _ (J  + 1)  (J  + 2) 


a(0  -*  2)  AE 


0 -♦  2 


2 (2  J + 1) 


( A — 3) 
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Expression  (A-3)  can  be  readily  shown  to  have  the  values  1.0,  1.0,  1.2, 
1.43  and  1.67  for  J = 0 - 4 respectively. 

The  conclusion  is  that,  at  temperatures  where  only  the  J = 0 and  1 
populations  are  significantly  populated,  the  degree  of  rotational  excitation 
will  be  independent  of  the  relative  populations  of  the  J = 0 and  1 states  for 
discharge  parameters  corresponding  to  characteristic  energies  well  above 
threshold.  This  observation  was  borne  out  by  detailed  discharge  predictions 
as  discussed  in  the  main  text.  This  behavior  may  not  hold  true  when  states 
with  J > 2 are  significantly  populated,  as  in  the  case  of  at  T = 300°K. 

In  this  limit  however,  the  effect  of  collisions  of  the  second  kind,  i.  e. 

J ■*  J - 2,  must  also  be  considered. 

At  present,  one  must  rely  on  theory  for  the  specification  of  the 
rotational  excitation  cross  sections  for  Ji  2,  For  purposes  of  complete- 
ness a set  of  qua  si-analytical  cross  sections  have  been  developed  for  use 
in  specifying  the  rotational  excitation  of  states  J = 0 - 3 in  H£  and  J = 0 - 4 
in  These  results  are  shown  in  Figs.  A-3  and  4.  The  solid  portions 

of  the  cross  sections  shown  are  as  determined  from  the  theoretical 

6 A 7 

predictions  of  Lane  and  Geltman  and  Henry  and  Lane.  ' The  dashed 
portions  represent  empirical  extrapolations  of  this  theoretical  work  which 
were  developed  in  the  present  analysis.  If  it  were  desired  to  examine  the 
effect  of  rotational  temperature  on  rotational  excitation  in  more  detail, 
such  cross  sections,  modified  slightly  to  improve  the  agreement  between 
theory  and  experiment,  could  be  used  in  conjunction  with  the  Boltzmann 
computer  code.  Such  an  effort  was  not  warranted  under  the  present 
program. 
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Electron  Excitation  Cross  Sections  for  Rotational  Transitions 
of  H~  as  a Function  of  Electron  Energy. 
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APPENDIX  B 


Rotational  Excitation  of  HC1 


For  most  molecules  a large  number  of  rotational  levels  can  be 
populated  at  room  temperature.  Because  of  this,  when  modeling  rotational 
excitation  by  electrons  one  must  consider  the  effect  of  collisions  of  the 
second  kind,  i.e.  electron-molecule  collisions  in  which  the  rotational 
state  of  the  molecule  is  decreased  rather  than  increased.  This  topic  has 
been  discussed  by  Hake  and  Phelps^  in  some  detail.  For  polar  molecules 
the  cross  section  for  rotational  excitation  is  given  by  9 


°J,  J + 1 (S) 


(J  + l)  Ry  a 

1 

(2  J+  l)  e 


- (e  -e 


^3 


(B  - 1 ) 


where  J is  the  rotationaj  level  of  the  molecule,  e is  the  electron  energy, 

Ry  is  the  Rydberg,  ar  - 8tth^  aQ^/3,  p is  the  electric  dipole  moment  in 
units  of  eaQ,  aQ  is  the  Bohr  radius  and  £j  is  the  energy  lost  by  the  electron 
in  the  collision,  2 (J  + 1)  B , where  BQ  is  the  rotational  constant  for  the 
molecule.  The  cross  section  for  rotational  de -excitation  may  be  deter- 
mined from  Eq.  (B-l)  by  the  use  of  detailed  balancing  and  is  given  by 


aj  + l,  J(S)  = 


(J  + l)  Ryar 

(2  J + 3)  e 


- . ,1/2*  1/2- 

(e  -t  €j)  + e 


[|s  + ej 


)1/2  -s1/2] 
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where  £j  is  now  energy  gained  by  the  electron  in  the  collision. 

To  properly  account  for  rotational  excitation  the  cross  sections 
(B-l)  and  (B-2)  should  be  included  for  each  rotational  level  J which  is 
populated  at  the  rotational  temperature  of  interest.  In  practice  this  is  quite 
unwieldy  and  effective  single  level  approximations  have  been  developed.  1 0 
The  resulting  effective  cross  section  is  given  by® 


'01 
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1 - exp  (-2B  /kT) 
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and  B is  an  effective  rotational  constant  equal  to  (B0kT)  . Shkarofsky 
et  al^  have  shown  that  the  average  energy  lost  in  a rotational  collision 
including  the  effects  of  rotational  de-excitation  is 


As  _ (A  E) 
e k T e 


(B-5) 


where  AE  is  2JB  . With  the  definition  of  B given  above  relationship  (B-5) 
reduces  to 


Ae  = 2 Bq  (B-6) 

Relationships  (B-3)  - (B-6)  were  used  in  the  present  analysis  in  the 

manner  prescribed  by  Shkarofsky  et  al.  10  which  is  slightly  different  than 

that  developed  by  Hake  and  Phelps.  8 Both  analyses  assume  that  2B'<sckT 

and  that  the  characteristic  energy  of  the  discharge  is  much  higher  than  the 

threshold  for  rotational  excitation.  In  the  case  of  HCl.which  has  an 

1 ft 

electric  dipole  moment  of  1.  08  x 10"  ° esu-cm  and  a rotational  constant  of 
1.  3 x 10-3  ev,  the  quantity  2B1  evaluated  at  300°K  has  a value  of  1.  2 x 10~2 
eV,  and,  thus,  the  above  approximations  are  only  marginally  valid. 

Finally,  the  computer  code  used  in  the  analysis  as  developed  by 
Yos  and  co-worker  si  1 had  assumed  that  the  logarithmic  factor  in  Eq.  (B-3) 
could  be  treated  as  a constant.  The  present  program  was  modified  to 
include  the  correct  functional  form  of  the  cross  section.  It  was  noted 
however,  that  the  discharge  predictions  performed  at  the  lowest  values 
of  E/N  considered,  did  not  accurately  satisfy  energy  conservation,  i.  e. 
the  balance  between  energy  transferred  to  electrons  from  the  electric 
field  and  energy  transferred  from  the  electrons  to  the  gas  molecules. 

It  is  suspected  that  this  inaccuracy  results  from  numerical  difficulties 
in  integrating  the  rotational  excitation  cross  section  over  the  electron  velocity 
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distribution  at  low  values  of  6K.  Inasmuch  as  the  predictions  at  low 
values  of  E/N  were  not  of  interest  for  laser  applications  this  matter 
was  not  pursued  further. 
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APPENDIX  C 

A Note  on  Vibrational  Excitation  of  HCL 
By  Electron  Impact* 


One  of  the  present  authors  recently  co-authored  a paper  describing 
an  experimental  determination  of  the  rate  constant  for  electron  impact 
vibrational  excitation  (EIVE)  of  HC1.  ^ In  this  experiment  an  electron 
beam  sustained  discharge  was  used  to  excite  gas  mixtures  of  HCI-CO-N2- 
Ar  and  the  EIVE  rate  constant  of  HC1  relative  to  that  for  CO  was  deter- 
mined via  infrared  fluorescence  measurements.  This  note  describes  a 
re-analysis  of  this  experimental  data  which  tends  to  reduce  the  uncertainty 
in  the  previously  published  rate  constants. 

12 

It  was  shown  earlier  that,  with  the  E-beam  sustainer  on,  the 
time  evoluation  of  vibrationally  excited  HC1  could  be  approximately 
represented  by  the  relationship 


[HC1*]  = lHCl*l 


1 - exp  ( — t / T ) ) 


where  [HC1  ] 


is  a quantity  proportional  to  the  EIVE  rate  constant  for 


HC1  and  T is  the  time  constant  for  collisional  quenching  of  HC1.  ' It  was 
further  shown  that  upon  termination  of  the  sustainer  pulse,  at  time  t = 
the  decay  of  HC1*  could  be  described  by  the  relationship 


[hci*1  = [hci*1  - ([HC1*] 


- |HC1>,=  1 ) exp  |(t  -t)/Tl 
td  ° 


(C-2) 


Authored  by  G.  E.  Caledonia  of  PSI  and  R.  E.  Center,  Mathematical 
Sciences  Northwest,  Inc.  , Bellevue,  Wash.  98009  and  published  in 
J.  Chem.  Phys.  64,  4237  (1976). 
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where  [HC1  |„  is  a residual  level  of  HC1  excitation  resulting  from  the 
vibration- vibration  exchange  reaction  with  N2,  and,  is  linearly  proportional 
to  the  EIVE  rate  constant  for  N2. 

Relative  measurements  of  the  fluorescence  histories  of  both  HC1 
and  CO  were  performed  in  this  experiment  for  a range  of  mixture  ratios 
and  discharge  parameter  E/N.  It  was  shown  that  the  EIVE  rate  constant 
for  HC1  could  be  determined  relative  to  that  for  CO  through  us"  of  the 
parameter  T (see  Eq.  13  of  Ref.  12)  which  could  be  deduced  directly  from 
the  HC1  fluorescence  measurements.  Unfortunately,  it  was  found  that  the 
observed  rise  time  of  the  HC1  emission  was  always  shorter  than  the  decay 
time,  in  some  instances  by  as  much  as  a factor  of  two.  Although  the 
measured  decay  time  was  more  in  accord  with  kinetic  estimates  for  T, 
both  measurements  were  used  to  provide  upper  and  lower  bounds  for  the 
rate  constant. 


It  has  since  been  realized  that  the  apparent  difference  in  the  rise 
and  decay  times  of  HC1  fluorescence  arose  from  an  error  in  the  data 
analysis.  In  the  evaluation  of  the  risetime  it  was  assumed  that  |HCl  '| 

D 

was  equal  to  [HC1  |s  s>  . This  was  not  the  case.  Indeed  for  the  range 

of  parameters  considered  [HC1*L  varied  from  50  to  90  percent  of 

t = T 

sjc  Xj 

[HC1  ls  s • The  effects  of  this  error  were  that  the  reduced  fluorescence 
data  appeared  to  be  qua  si- exponential  and  that  the  deduced  risetime  was 
always  less  than  the  actual'value  of  T.  On  a proper  evaluation  of  the  data 
employing  a detailed  kinetic  analysis,  it  was  found  that  TriSe  was  indeed 
equal  to  Tdecay.  Thus  the  fluorescence  data  may  be  consistently  explained 
by  the  lower  value  of  the  rate  constants  shown  in  Fig.  3 of  Ref.  12. 


Furthermore,  it  has  been  found  that  the  EIVE  rate  constant  may 

also  be  determined  from  the  data  without  reference  to  the  CO  fluorescence. 

It  may  be  shown  that  the  ratio  of  the  EIVE  rate  constant  for  HC1  relative 

to  that  for  N2  can  be  determined  from  the  measured  ratio  of  |HC1  | t 

$ ^ 
to  [HC1  ]„  (see  Eqs.  11  and  14  of  Ref.  12).  The  rate  constants  deduced 

in  this  manner  were  typically  20%  higher  than  those  deduced  from 

comparison  with  the  CO  fluorescence  measurements.  This  additional 

comparison  provides  a good  consistency  check  on  the  experiment. 

The  quantity  measured  in  this  experiment  is  the  EIVE  rate  constant 
for  HC1  vs.  the  characteristic  energy  of  the  discharge.  This  measured  rate 
constant  is  thus  a function  of  the  specific  non-Maxwellian  electric  energy 
distribution  peculiar  to  the  experimental  range  of  E/N  and  mixture  ratio. 
Indeed  the  rate  constants  shown  are  not  necessarily  representative  of  those 
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for  other  gas  mixtures  under  discharge  conditions  corresponding  to  the 
same  characteristic  energy.  Such  information  can  only  be  determined 
from  the  energy  dependent  EIVE  cross  section  for  HC1  which  cannot  be 
deduced  from  the  present  data. 

There  have  been  two  recent  measurements13,  14  of  the  EIVE  cross 
section  for  HC1.  In  the  first  of  these13  the  energy  dependence  of  the  cross 
section  was  measured  in  the  energy  range  of  0 - 10  eV.  This  was  a relative 
measurement  and  was  normalized  at  an  electron  energy  of  10  eV  to  an 
absolute  measurement  of  the  total  scattering  cross  section  for  HC1  taken 
by  Bruche.  The  peak  cross  section  observed  was  1.  3 x 10"1®  cm2  ± 50% 
where  the  error  bars  represent  the  experimental  uncertainty  but  do  not 
include  the  uncertainty  in  Bruche's  cross  section. 

The  EIVE  cross  section  for  HC1  was  also  studied  by  Ziesel  et  al14 
in  the  onset  region.  These  measurements,  taken  between  0.  36  - 0.  5 eV 
were  consistent  with  the  observations  of  Rohr  and  Linder;  however,  once 
again  the  measurements  were  relative,  in  this  instance  normalized  to  the 
dissociative  attachment  cross  section  for  HC1.  The  peak  cross  section  was 
estimated  to  be  in  the  range  of  2 x 10- - 2 x 1 0*  1 5 cm2  depending  upon 
the  choice  of  the  dissociative  attachment  cross  section.  H 

Thus  there  is  still  a significant  uncertainty  in  the  magnitude  of  the 
EIVE  cross  section  for  HC1.  Discharge  predictions  have  been  performed 
for  the  experimental  conditions  of  Ref.  12,  using  the  relative  HC1  cross 
section  determined  by  Rohr  and  Linder13  in  order  to  ascertain  what 
magnitude  cross  section  would  be  consistent  with  the  experimental  rate 
constants.  It  was  found  that  a peak  cross  section  of  0.  8 x 10-13  cm^ 
provided  good  agreement  with  the  measured  rate  constants.  This  value 
is  within  the  error  bars  of  the  other  measurements13,  14  and  thus  it 
would  appear  that  the  peak  cross  section  for  EIVE  of  HC1  is  approximately 

1 x 10"15  cm2.  r y 

There  is  one  last  point  to  be  made  concerning  the  results  presented 
in  Ref.  12.  In  that  work  it  was  found  that  approximately  20%  of  the  observed 
HC1  fluorescance  arose  from  transitions  higher  than  v = 1 -*  0 and  it  was 
suggested  that  this  emission  resulted  from  direct  electron  excitation  into 
the  second  or  higher  vibrational  levels.  *2  Since  that  time  detailed 
anharmonic  oscillator  kinetic  calculations,  including  all  relevant  inter- 
mode and  intra-mode  vibration- vibration  reactions  for  HC1  and  N2,  ^ 
have  been  performed  for  typical  experimental  discharge  and  mixture 
conditions.  These  calculations  demonstrate  that  approximately  10%  of 
the  vibrationally  excited  HC1  could  reside  in  states  higher  than  v = 1 even 
if  there  were  no  direct  electron  excitation  of  these  states.  This  result 
is  solely  due  to  vibration- vibration  exchange  collisions  occurring  over  the 
60n  sec  duration  of  the  discharge.  Thus  no  information  concerning  EIVE  to 
vibrational  levels  greater  than  v = 1 can  be  deduced  from  the  data. 
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In  the  second  experiment,  HCl/inert  gas  mixtures  are  directly 
excited  by  a pulsed,  high  current  e-beam  sustainer  discharge.  The 
discharge  conditions  were  modeled  by  developing  a set  of  electron  impact 
cross  sections  for  HCl.by  a critical  review  of  existing  data.  ^The  results 
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